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WITH  APPLICATION  TO  SOME  RE-ENTRY  EXPERIMENTS 

by 

A.  P.  Waterfall 


SUMMARY 

An  automatic  technique,  utilising  a  large  high  speed  computer  for  the 
analysis  of  records  obtained  from  accelerometers,  pressure  gauges  and  gyros 
installed  in  the  re-entry  vehicles  of  the  Black  Knight  rocket,  is  described. 

It  is  shown  that  this  enables  a  suitable  arrangement  of  only  three  accelerometers 
to  yield  trajectory  data,  angular  motion  and  the  aerodynamic  characteristics  of 
the  re-entry  body.  Examples  of  the  analysis  of  accelerometer  and  rate  gyro  data 
are  given  and  the  lines  of  further  development  and  some  applications  are 
indicated. 
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1  INTRODUCTION 

Black  Knight  is  a  ballistic  research  test  vehicle  primarily  designed  for 
atmospheric  re-entry  studies  and  firings  at  Woomera  have  continued  since  1958* 
An  important  part  of  the  Black  Knight  programme  is  the  on-board  measurement  of 
the  dynamic  motion  and  aerodynamic  characteristics  of  the  re-entiy  vehicles 
flown.  These  re-entty  vehicles  are  axially  symmetric  and  mainly  oonical  but 
some  are  spherical  in  shape  and  re-enter  the  atmosphere  at  speeds  of  up  to 
16000  ft/sec.  The  current  Black  Knight  programme,  Project  Dazzle  is  a  re-entiy 
physics  programme  measuring  the  radar  and  optical  characteristics  of  re-entry 
vehicles.  Measurements  by  ground  based  instrumentation  of  either  optical  or 
IR  radiation  or  radar  returns,  due  to  wake  or  ionisation  effects  associated 
with  the  re-ontry  body,  are  strongly  affeoted  by  body  incidence.  In  these 
experiments  it  is  therefore  important  to  correlate  the  angular  motion  of  the 
vehicle  with  the  ground  measurement,  as  well  as  continuing  the  basio  measure¬ 
ment  of  aerodynamic  characteristics.  These  more  exaoting  requirements  plus  a 
higher  firing  rate  made  it  desirable  to  deviBe  some  form  of  automatic  analysis 
scheme  for  the  re-entry  dynamics  data,  especially  for  that  obtained  from 
instruments  carried  by  the  re-entiy  vehicle.  The  techniques  of  measurement 
and  data  analysis  that  are  under  development  form  the  substance  of  this  paper. 

There  are  two  possible  approaches  to  making  aerodynamic  measurements  in 
free  flight.  The  first  is  essentially  similar  to  wind  tunnel  techniques  in 
which  forces  and  moments  are  measured  directly  along  with  incidence  a,  so  that 
force  and  moment  coefficients  C^,  and  may  be  calculated  as  functions  of 
a.  The  alternative  approach  arrives  at  values  of  the  aerodynamio  derivatives 
CNa’  CMq  by  studying  the  tynamic  response.  This  does  have  limitations  in 
the  general  case,  because  the  form  of  the  coefficient  vs.  incidence  curve  cannot 
be  completely  determined.  However,  adequate  results  are  usually  obtained  by 
assuming  some  polynomial  form  to  the  curve.  The  advantages  of  the  dynamic 
method  is  that  it  needs  much  less  instrumentation  than  the  direct  method,  and 
it  is  much  in  favour  where  space  is  at  a  premium.  It  is  this  type  of  method 
with  which  we  are  primarily  concerned. 

Analysis  of  dynamic  response  always  needs  a  theoretical  dynamic  model. 

The  one  used  for  the  analysis  of  the  ’Black  Knight'  data  is  discussed  in 
Section  2,  together  with  an  outline  of  the  theory  of  re-entry  and  an  appraiaal 
of  some  of  its  limitations.  This  leads  to  some  theoretical  considerations, 
under  Section  3,  of  the  kind  of  instrumentation  that  may  be  carried  by  re-entry 
vehicles.  Accelerometers  in  particular  are  considered  in  detail,  because  they 
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are  usually  the  simplest  and  most  compact  system.  It  is  shown  how  some 
inherent  difficulties  may  be  overcome  and  how  three  accelerometers,  suitably 
deployed,  can  monitor  almost  all  the  data  necessary  for  a  complete  dynamic 
analysis. 

Section  3  also  gives  a  brief  account  of  the  'Black  Knight'  Re-entry 
Vehicle  instrumentation  system.  In  the  most  comprehensive  system,  in  which 
direct  measurements  of  coefficients  are  also  attempted,  rate  gyros, 
differential  pressure  systems,  accelerometers  and  ground  based  ballistic 
cameras  yield  the  basic  dynamics  data.  However,  the  usual  practice  is  to 
limit  the  on-board  instrumentation  to  accelerometers  only,  because  of  shortage 
of  space. 

There  are  two  methods  of  dynamic  analysis  of  re-entry  data:  tho  first 
is  suitable  for  a  graphical  analysis  by  hand,  the  second  for  an  automatic 
method  on  a  computer.  The  graphical  method,  described  in  detail  under 
Appendix  A  is  essentially  one  of  progressive  simplification  of  the  data. 

It  starts  with  an  estimate  of  the  trajectory  30  as  to  obtain  dynamic  pressure 
as  a  function  of  time.  The  various  components  of  the  dynamic  motion  are 
then  separated  and  progressively  reduced  to  a  plot  of  the  motion  and  estimates 
of  the  aeroctynamic  derivatives. 

It  may  be  possible  to  adopt  the  graphical  method  for  computer  use  but 
it  would  be  most  untidy;  it  is  better  to  adopt  a  different  approach  based 
on  least  squares  fitting  methods.  Hie  basic  idea  is  to  begin  by  making  a 
guess  at  the  parameters  in  the  dynamic  model,  which  is  then  used  to  produce 
theoretical  numbers  for  comparison  with  the  actual  data.  The  method  of 
differential  corrections,  described  in  Section  4,  is  then  used  to  improve 
the  estimates  of  the  parameters  until  the  best  fit  to  the  data,  on  the  least 
squares  criterion,  is  obtained.  Such  a  technique  is  ideally  suited  to  com¬ 
puter  use,  and  it  has  the  great  advantage  that  the  model  may  be  as  compre¬ 
hensive  as  we  please,  whereas  there  are  practical  limits  to  a  graphical 
analysis  performed  by  hand. 

Automatic  methods  using  a  digital  computer  were  developed  first  for 
the  analysis  of  trajectory  data  from  ballistic  cameras  and  accelerometers. 
These  are  briefly  described  under  Section  5.1  and  illustrated  with  actual 
results  in  which  a  smoothed  trajectory  is  obtained  along  with  estimates  of 
as  a  funotion  of  height.  Such  least  squares  methods  have  long  been  used 
for  analysing  the  dynamio  motion  in  ballistic  ranges  to  obtain  the  aerodynamic 
derivatives.  The  dynamic  model  in  this  case  is  fairly  simple.  The  model 
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described  in  Section  2  is  necessarily  more  complex  to  take  account  of  the  large 
changos  in  environment  during  re-entry  and  a  large  high-speed  computer  now 
becomes  essential  if  the  computing  time  is  to  be  kept  to  a  reasonable  duration. 
These  programmes,  identical  in  principle  but  much  more  complex  than  the  tra¬ 
jectory  programmes  are  described  briefly  under  Section  5*2.  typical  results  w® 
also  discussed  showing  how  the  angular  motion  is  fitted  and  estimates  made  of 
aerodynamic  derivatives. 

Operating  details  of  the  ano’/sis  programmes  are  given  in  Appendix  C. 

The  report  ends  with  a  discussion  of  possible  future  developments  and 
applications. 

2  THE  DYNAMIC  MODEL 

It  is  convenient,  for  the  purposes  of  analysis,  to  consider  re-entry 
dynamics  in  two  parts,  and  to  take  first  the  motion  of  a  point  having  given 
aerodynamic  properti; s ,  and  then  the  motion  of  the  body  about  its  C  of  G  as  it 
descends  through  the  ahnosphere.  The  first  loads  to  a  description  of  the 
trajectory  of  the  body  in  terms  of  height  and  velocity  as  a  function  of  time, 
and  the  latter  to  the  dynamic  motion.  Strictly  speaking  the  trajectory  and 
motion  are  inter-related,  but  the  technique  described  in  this  paper  is  concerned 
with  re-entry  experiments  in  which  the  angles  of  incidence  are  intended  to  be 
small.  In  such  cases  the  dynamic  motion  produces  only  small  perturbations  of 
the  trajectory  which  can  usually  be  determined  quite  accurately  by  making  only 
very  rough  allowances  for  the  effect  of  the  lateral  motion. 

2.1  Equations  of  motion  of  a  particle  over  a  fixed  earth  in  two  dimensions 

In  the  equations  of  motion,  effects  arising  from  the  rotation  of  the 
earth  are  ignored  but  gravity  terms  are  included  although  even  these  are 
negligible  in  comparison  with  the  aerodynamic  terms  except  during  the  very 
early  stages  of  re-entry.  The  equations  of  motion  then  take  the  following 
convenient  form  for  small  values  of  incidence: 


e 

V  = 

(-  D  -  Na)/m  -  g  sin  0 

(D 

ve  = 

(-  Da  +  N)/m  -  g  cos  0  +  cos  0/r 

(2) 

II 

•A 

r  =  V  sin  0  • 

(3) 

2 

Where  all  parameters  axe  defined  in  Fig.1  and  where  N  a  ?pV  sci  C„, 

2  2  2 

D  s  ?  pV  s  and  g  =  gQ  ryr  and  p  is  air  density,  s  the  reference  area, 

C.  axial  force  coefficient  and  C„  the  normal  foroe  derivative.  Note  that  D 
A  Na 

rather  than  A  is  used  for  axial  force  in  this  paper. 
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By  making  certain  assumptions  about  the  variation  of  p  with  height, 
ant  ignoring  all  except  the  drag  term,  the  above  equations  may  bo  integrated 
to  obtain  expressions  for  velocity  and  height  as  a  function  of  time.  However, 
the  advent  of  very  high  speed  computers  makes  these  rather  inflexible 
trajectory  models  unnecessary.  Much  more  useful  results  may  be  realised 
by  integrating  the  equations  by  a  numerical  method.  It  is  then  possible, 
for  example,  to  determine  the  variation  of  with  Mach  number  and  altitude 
rather  than  the  mean  ballistic  coefficient.  This  is  discussed  further  in 
Sections  3  and  1*., 


2.2  The  three-dimensional  gyrations  of  a  missile  descending  through 
the  atmosphere 

The  basic  equations  of  motion  of  a  spinning  body  during  re-entry  have 
much  in  common  with  the  stability  equations  of  aircraft  and  spinning  shells, 
with  the  extra  complication  of  rapidly  varying  dynamic  pressure  and  aero¬ 
dynamic  coefficients.  Theoretical  studies  of  the  dynamics  of  re-entry  have 
received  a  great  deal  of  attention  of  recent  years,  although  many  of  these 
have  made  simplifying  assumptions  which  render  them  unsuitable  for  our 
purpose.  The  approach  given  in  Ref.1  is,  perhaps,  tho  most  comprehensive, 
and  has  been  followed  in  tho  authors  work.  Only  relevant  results  will  bo 
quoted  and  these  without  proof. 


We  are  concerned  with  small  angular  perturbations  of  a  missile  about 
the  trajectory.  Small  perturbations  are  a  necessary  assumption  of  the 
mathematical  solution,  but  it  is  not  a  serious  restriction  as  incidences  of 
up  to  40  degrees  are  allowed.  It  is  assumed  that  the  basic  unperturbed 
trajectory  is  a  straight  line.  The  effect  of  perturbations  is  to  disturb 
both  the  trajectory  and  the  orientation  of  the  missile  axis.  To  specify 
these  angular  perturbations  we  will  use  an  orthogonal  reference  frame  fixed 
in  the  body.  (N.B.  thi3  differs  from  Ref.1  whore  the  axes  are  fixed  in 
space.)  We  take  X  along  the  axis  of  symmetry  of  the  missile  and  the  origin  0 
at  the  centre  of  gravity.  Then,  looking  in  the  direction  of  motion  from 
behind  the  C  of  G-,  the  angular  perturbations  can  be  represented  as  projec¬ 
tions  on  the  YZ  plane  thus:- 
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where  the  missile  has  perturbations,  a  (pitch),  p  (yaw)  and  angular  velocity 
resolutes  P  (spin),  q  (pitch),  r  (yaw). 


The  equations  of  motion  in  terms  of  the  complex  quantities  a -ip  are 
derived  in  Ref.1.  Certain  terms  are  negligible  from  a  practical  point  of 
view.  On  deleting  such  terms,  the  following  equation  remains:- 

^  *  3^[ipv,  -  iia)  ♦  air  (i  -  4)] 


o  -  ip  r~ 


i  P  v2sec, 


Ma  _  p2 


H) 


i  P  v2 ,  t  c*a 


a  -  ip 

e  e 


(0 


This  is  the  body  axes  counterpart  of  equation  (12)  of  Ref.1,  with  the  addition 
of  an  aerodynamic  asymmetry  term  on  the  right  hand  side.  It  should  bo  noted 
that  a  -  ip  is  small  compared  with  a  -  ip  so  that  terms  involving  V  have  been 

©  ©  p 

omitted  from  the  right  hand  side  in  comparison  with  v  term. 

This  equation  (4)  could  be  used  as  it  stands  as  the  dynamical  model  for 
the  dynamic  motion  in  the  same  way  as  we  do  for  the  trajectory.  This  would  bo 
a  retrograde  step,  however,  as  a  great  deal  of  pl\ysical  insight  would  be 
lost  -  the  equation  is  much  too  complex,  which  is  not  true  in  the  case  of  the 
simple  trajectory  equations.  An  analytical  solution  to  equation  (4)  is  there¬ 
fore  required. 

The  solution  may  be  considered  as  the  sum  of  two  particular  integrals, 
that  of  motion  about  trim  and  that  of  trim  itself.  !3he  first  is  obtained  by 
making  aQ  -  ipg  zero  in  free  space,  a -ip  taking  the  free  space  value.  Ihe 
second  is  obtained  by  making  a -ip  and  its  derivatives  zero  in  free  space  and 
studying  the  effect  of  a  -  ip  as  the  body  enters  the  atmosphere.  The  latter 
is  a  complex  problem  and  will  be  referred  to  later.  The  solution  offered  here 
is  for  the  case  of  zero  trim  or  that  part  of  ro-entiy  well  away  from  the 
resonance  region.  A  quasi-static  solution  to  equation  (4)  for  the  motion  about 
trim  is  derived  in  Rof.1.  In  terms  of  body  axes  it  is  as  follows: 
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a  -10 


whero 


«  +  AP/2B  -  P) 


+  S  oxp 


t 

i  I  (uJ  +  AP/2D  -  p) 
0 


•  *  • 


(5) 


w2  =  -  7  p  V2  s  C  CUa/B  +  A2P2/4B2 


R 


S 


a 


b 


o 


![cl/-*2  V6] 

Vs] 


and  u)  ,  R  and  S  represent  the  situation  at  t  =  0. 

0  0  o 

This  is  the  well  known  form  of  two  arms  rotating  in  opposite  directions, 
the  R  arm  against  the  spin  at  a  rate  (-  W  +  AP/2B  -P)  and  the  S  arm  with  the 
spin  at  a  rate  (w  +  AP/2B-P).  In  the  absence  of  spin  the  arms  rotate  at 
equal  and  opposite  rates,  and  so  would  draw  an  ellipse  at  rate  u>.  The  effect 
of  spin  is  to  cause  this  ellipse  to  precess  at  rate  P(l-A/2B)»  resulting  in 
the  familiar  flower  petal  pattern  of  Fig. 4. 

The  limitations  of  this  model  are  fully  discussed  in  Ref.1  and  it  would 
appear  that,  for  all  practical  purposes,  it  is  valid  for  symmetrical  missiles 
with  linear  aerodynamic  properties  below  200  000  ft.  Above  that  height  it 
would  be  seriously  in  error  if  the  spin  were  very  small,  much  smaller  than 
would  be  desirable  in  practice.  The  two-arm  pattern  also  romains  qualitatively 
true  in  the  presence  of  non-linear  aorodynamios  and  the  sort  of  small  asym¬ 
metries  which  are  the  result  of  manufacturing  imperfections.  Allowance  for 
some  of  these  effects  will  be  discussed  in  Section  6.  Howevor,  experience 
has  shown  that  this  linear  symmetric  theory  provides  a  very  useful  dynamical 
model  for  a  preliminary  analysis  of  re-entry  motion. 
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Before  considering  methods  of  measurement  it  will  be  helpful  to  write  down 


the  derivatives  of  a -ip.  Define: 

t 


Yr  =  J  ("  w  +  AP/2B  -  P)  dt  + 


Yr 


(6) 


Yg  =  j  (5  +  AP/2B  -  P)  dt  +  Yg 


(7) 


then  from  equation  (5) 


Now  in  body  axes 


and 


substituting  we  have 


i  i 

„  YS 


a  -ip  =  Re  +  S  e 


q  +  ir  -  -  a-  ip-iPa-ip 


q  +  ir  =  -a-ip-iPa-ip 


(8) 

(9) 

(10) 


q+ir  =  -  i[-  »  +  AP/2B]  Re  K  -  i[»  +  AP/2B)  S  e 

iy 

— ~  =  [-  »  +  AP/2B][-  «  +  AP/2B  -  P]  R  e  R 


YS 


+  [»  +  AP/2B][»  +  AP/2B  -  P]  S  e 


writing  out  in  full,  and  including  constant  trim  angles  a.  and  p  we  obtain 

t  V 


a  =  R  cos  Yu  +  8  cos  Ye  +  a, 

K  5  X 

P  =  -  R  sin  yr  -  S  sin  Yg  +  Pt 


] 


(11) 


q  =  R  [-  »  +  AP/2B]  sinYR  +  S[£>  +  AP/2B]  sinys 
r  b  -  R  [-  u  +  AP/2B]  cosyR  -  S[u>  +  AP/2B]  cos  y® 


} 


(12) 
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4  =  R  [-tt>  +  AP/2B][-Jj  +  AP/2B-P]  003  yR  +  S[w  + AP/2B][tt)  + AP/2B-P]  cosy  S 

r  =  R  [-  u>  +  AP/ 2B] [-  io  +  AP/2B  -  P]  sin  yR  +  S[»  +  AP/2B][u>  +  AP/2B -P]  sin  y« 

...  (13) 


3  METHODS  OP  MEASUREMENT 


Instrumentation  for  the  monitoring  of  a  particular  flight  is  divided 
between  on  board  instruments  and  ground  based  cameras  or  radars.  Hie  latter 
are  required  to  monitor  position  with  time  and  will  not  be  discussed. 
Trajectory  analysis  normally  consists  in  combining  ground  data  and  airborne 
accelerometer  data  to  arrive  at  a  best  fit  to  the  data.  Accelerometers, 
however,  are  affected  by  the  oscillations  of  the  body  and  need  correcting 
before  use  in  a  trajectory  analysis.  It  will  therefore  be  appropriate  to 
consider  methods  of  monitoring  the  dynamic  motion  with  on  board  instruments 
before  the  monitoring  of  the  trajectory. 

Three  instruments  are  in  common  use  for  the  measurement  of  dynamic 
behaviour,  acoelerometers,  rate  gyros  and  differential  pressure  devioes. 

Rate  gyros  monitor  q+ir  directly  and  are  perhaps  the  most  useful  and  accurate 
instrument  we  can  use.  Differential  pressure  measurements  give  a—i{3 
assuming  dynamic  pressure  is  known.  They  need  pre-flight  calibrations  and 
sometimes  suffer  from  slow  response.  Both  the  latter  are  well  understood 
and  will  not  be  discussed  further.  Acoelerometers,  on  the  other  hand,  are 
used  more  than  any  other  instrument  in  ballistic  missiles  because  of  their 
small  size  and  modest  power  requirements.  Unfortunately  they  are  also  the 
least  understood,  and  give  the  most  opportunity  of  being  incorrectly 
interpreted.  They  are  therefore  worthy  of  careful  consideration. 

3«1  The  use  of  accelerometers 


The  usual  purpose  of  accelerometers  is  to  measure  the  acceleration  of 
the  centre  of  gravity,  and  for  this  they  must  be  located  at  the  centre  of 
gravity.  The  acceleration  along  the  three  principal  axes  (A^,  A  ,  A^)  is 
given  by: 
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Ax  =  “  ipV2  a  C^/m 

Ay  * 

\  =  ipV2  s  C„a(a-<it)/m  . 


A  is  the  most  used  in  trajectory  determination.  A  and  A  can  be  used  in  the 
x  y  2 

same  way  as  differential  pressure  measurements  to  determine  a  and  p,  or,  if 

these  are  known  to  find  C„  . 

Net 

As  it  is  frequently  impossible  to  locate  accelerometers  at  or  even  near 
the  C  of  G,  we  shall  now  consider  the  acceleration  measured  at  a  general 
point.  As  before,  let  the  principal  axes  be  CXYZ  where  0  is  at  the  C  of  &. 

Let  the  resolutes  of  velocity  be  (u,  v,  w)  and  let  the  resolutes  of  angular 
velocity  Q  be  (P,  q,  r).  We  consider  a  general  point  L  located  at  (x,  y,  s). 


The  velocity  of  L  in  body  axes  is  given  by: 

V  =  [»»*)♦  p  $  l] 

=  [u+qz-ry,  v-Pz  +  rx,  w+iy-  qx]  . 
The  acceleration  of  L  is  ^  +  [Hv]  which  is 

=  [u  +  qz  -  ry,  v  -  Pz  +  fx,  w  +  Fy  -  qx]  ♦ 


r  P  q  r  1 

|_u  +  qz  -  ry  v-Pz  +  rx  w+]ty  -  qxj 

=  [(u  +  qw  -  rv)  -  (q2  +  r2)  x  +  (qP  -  r )  y  +  (pr  +  q)  z, 

(v  +  ru  -  Pw)  +  (Pq  +  r)  x  -  (P2  +  r2)  y  +  (qr  -  P)  z, 

(w  +  Pv  -  uq)  +  (Pr  -  q)  x  +  (qr  ♦  P)  y  -  (q2  +  P2)  z]  . 
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The  first  terms  in  each  component  are  actually  the  C  of  S  acceleration 
in  body  axes  (Ax,  Ay,  Az) ,  so  that: 

Acceleration  parallel  to  X  axis 

ax  =  "  ip  V2  8  CA^m  “  ^q2  +  r2)  x  +  (p1  ”  **)  y  +  (r  +  q)  z  . 

...  (14) 

Acceleration  parallel  to  Y  axis 

ay  =  ipV2sCNa(p  -  0t)/m  +  (Pq  +  f)  x  -  (P2  +  r2)  Y  +  qrz  . 

...  05) 

Acceleration  parallel  to  Z  axis 

az  =  2  pv  s  CNa(a  -  at)/m  +  (Pr  -  q)  x  +  q  ry  -  (P  +  q  )  z 

...  (16) 

where  P,  q,  r,  a,  p,  q,  r  axe  given  by  equations  (ll)  to  (13)  and  £>  is 
assumed  to  be  negligible. 

It  will  readily  be  appreciated  that  the  general  situation  is  very  com¬ 
plicated  and  that  little  can  be  done  to  simplify  the  above  equations. 

2  2 

Fortunately  those  terms  involving  qr  and  q  and  r  ,  which  destroy  the  simple 
two  component  form  of  the  C  of  G  accelerations,  are  often  small  and  may  be 
neglected  in  analysis  by  hand.  Corrections  are  readily  made  for  centrifugal 
effects  in  automatic  analysis  techniques  when  accelerometers  have  been 
scattered  about  a  re-entry  vehicle  in  an  arbitrary  fashion,  but  are  a  nuisance 
particularly  in  the  case  of  trajectory  analysis.  However,  accelerometers  can 
be  made  to  yield  extra  useful  data  if  displaced  from  the  C  of  G  with  car®  and 
thought,  especially  axial  displacement  of  lateral  accelerometers.  This 
special  case  will  now  be  considered. 

3*2  The  use  of  axially  displaced  lateral  accelerometers 

Suppose  we  have  an  accelerometer  measuring  acceleration  parallel  to 
the  Y  axis  but  displaced  from  the  C  of  G  along  the  X  axis.  Wien 

ay  a  ip  V2®  CNa(p- Pt)/m  +  (Pq  +  r)  x  (17) 

substituting  for  p,  P,  q  and  r  from  equations  (ll),  (12)  and  (13) 
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fly  =  C-  5P  V2s  CNa/m  +  x(u>  -  AP/2B)2]  R  sin  Yr 

+  [-  ipV2s  CNa/m  +  x(u>  -  AP/2B)2]  S  sin  Yg  •  (18) 

Hence  the  general  effect  of  an  axial  displacement  is  to  change  tire  magnitude 
of  the  two  oscillatory  components  of  the  accelerometer  reading.  A  backward 
displacement  (-x)  increases,  a  forward  displacement  decreases  the  magnitude  of 
the  components. 

If  we  have  two  accelerometers  at  axial  stations  x1  and  x^,  and  subtract 
their  outputs,  we  eliminate  terms  involving  Cjj  and  are  left  with  the  centri¬ 
fugal  terms  only: 


A  =  (x1  -x2)  [ (<u-AP/2B)2  R  sin  Yr  +  (w  +  AP/2B)2  S  sin  Yg]  •  0?) 


From  a  plot  of  A  it  is  possible,  in  principle,  to  obtain  R  and  S  without  even 
a  knowledge  of  the  trajectory.  This  is  the  basis  of  the  method  suggested  by 
Nelson.  In  his  method,  an  orthogonal  pair  of  accelerometers  is  requirod  at 
each  axial  station  because  the  data  is  analysed  by  the  graphical  method  des¬ 
cribed  in  Section  4.1.  However,  the  automatic  method  of  the  author  can  fit  a 
curve  to  a  plot  of  A,  or  single  record  and  deduce  R  and  S  as  well  as  other 
parameters.  Moreover  can  be  deduced  by  substitution  of  an  analysis  of  (19) 
back  into  (l8).  Consequently,  with  a  third  accelerometer  to  measure  longi¬ 
tudinal  acceleration,  a  complete  dynamic  analysis  becomes  possible.  For,  is 
determined  from  the  trajectory,  from  u>,  C^  and  from  the  damping  of 

R  and  S.  In  addition  CT  is  also  given  by  the  difference  of  C„  and  C..  Thus 

La  Na  A 

three  accelerometers  can  be  deployed  so  that  all  the  aerodynamic  properties  of 
the  re-entry  vehicle  are  monitored.  The  axial  displacement  (x^  -  xR)  has  to  be 
large  and  the  accelerometers  need  to  be  very  accurate  to  give  superior  results 
to  a  combination  of  rate  gyros  and  accelerometers. 


It  is  worth  noting  that  over  most  of  the  re-entry  u>  is  very  much  bigger 
them  AP/2B  so  that  equation  (l8)  may  bo  reduced  to  the  form: 


a 

y 


v2  j 


(P-  Pt)0  x 


ch</b 


(20) 


The  term  at  x  C^/B  is  a  constant  (approximately),  so  that  the  effect  of 
moving  an  accelerometer  axially  from  the  C  of  &  is  to  change  the  magnitude  of 
oscillations  by  this  factor. 
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3.3  Determination  of  the  trajectory  from  accelerometers 

If  a  trajectory  is  to  be  integrated  it  is  necessary  to  determine  the 
aerodynamic  terms  of  equations  (l)  and  (2),  these  are  not  given  directly  by 
acceleration  measured  in  the  body.  The  term  (-  D  a+N)  of  equation  (2)  is 
oscillatory  and  it  has  been  found  reasonable  to  assume  that  its  net  effect 
on  0  is  zero.  The  problem  is  reduced  therefore  to  a  determination  of  the 
value  of  (D+Na)  of  equation  (l). 

The  major  part  of  this  term  is  measured  by  the  longitudinal  accelero¬ 
meter  but,  as  will  be  seen  from  equation  (llf)  there  are  complications  caused 
by  displacements  from  the  C  of  G.  It  may  be  shown  that  lateral  displacements 
result  in  symmetrical  oscillatory  errors,  which  may  be  removed  by  reading 
the  mean  value  of  acceleration.  This  leaves  the  effect  of  axial  displacement 
only  vis: 

ax  =  “  i  P  V2  s  CA/m  -  (q2  +  r2)  x  . 

Substituting  from  equation  (lO),  assuming  that  uS  >>  AP/2B,  and  taking  the 
mean  value  of  a^  we  have 

ax  =  -  i  pV2  sCA/m  -  u>2(R2+  S2)  x  . 

I" 2  2  - 

Here  I  R  +  S  is  the  rms  of  incidence  a  if  trim  is  assumed  to  be  zero. 
Substituting  for  u>  also  we  have 

ax  =  -  (C  -  mt  x  a2  C^B)  £pV2  s/m  .  (21 ) 

Combining  with  equation  (l)  we  have: 

V  b  a  -(l+mdx  C„  /B  C„  )  —  -  g  sin  0  .  (22) 

x  Mcc  Na  m 

Using  the  lateral  accelerometers  an  estimate  of  a  can  be  made  by  estimating 
the  rms  value  a  which  is  equal  to  N/m.  Thus: 

y 

V  =  ax  -  [1  4  me  x  CMo/B  C^]  a^p^sC^  -  g  sin  0  (23) 

and 

V0  b  ”  gQ  cos  6  +  V2  cos  0/r  (24) 

h  >  r  b  V  sin  0  .  (25) 
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Given  an  initial  height,  velocity  and  climb  angle  these  equations  may  be 

integrated  to  give  a  re-entry  trajectory,  using  values  of  a  and  a  from  the 

x  y 

accelerometer  records.  It  will  be  shown  in  Section  5*1  how  this  may  be  com¬ 
bined  with  a  least  squares  fit  to  ballistic  camera  observations.  Although  only 
approximate  allowance  has  been  made  for  the  effect  of  dynamic  motion  on  tho 
accelerometer  and  the  trajectory,  it  has  been  found  adequate  in  practice. 

3.4  Black  Knight  instrumentation  system 

A  considerable  quantity  of  sub-miniature  instrumentation  can  be  carried 
in  Black  Knight  re-entry  vehicles.  Typical  instrumentation  payloads,  i.e. 
excluding  the  outer  shell,  internal  support  structure  and  ballast  weights,  have 
been  up  to  60  lb  in  weight.  Output  from  transducers  can  be  transmitted  to 
ground  receivers  via  radio  telemetry  or  alternatively  recorded  on  magnetic 
tape.  Provided  the  recorded  tape  is  stored  in  an  armoured  cassette,  recovery 
presents  little  difficulty. 

The  latter  system  which  is  not  troubled  by  ionisation  blackout  during 
re-entry,  has  been  used  exclusively  for  a  number  of  years  and  a  range  of 
magnetic  tape  recorders  has  been  developed.  The  most  favoured  version,  the 
sub-miniature  recorder,  is  fully  armoured  and  provides  a  total  record  time  of 
3?  minutes  in  8  tracks  on  ?  inch  tape  running  at  1 7/8  inches  per  sec.  Hie 
tape  recorder  is  switched  on  shortly  before  re-entry  commences. 

One  tape  track  is  used  for  recording  time  in  the  form  of  pulses  at  a 
prf  of  2*5  kc/s  provided  by  a  crystal  controlled  source.  The  main  system  of 
recording  on  other  tracks  is  by  frequency  modulation.  The  input  to  the  modu¬ 
lator  is  ±2*5V  D.C.,  the  output  is  1*5  kc/s  ±600  c/s  feeding  direct  to  the 
recording  heads.  Each  F.M.  input  is  time  multiplexed  by  a  24  way  switch 
running  at  3  revs  per  sec.  The  principle  use  of  multiplexing  on  channels 
recording  dynamics  information  is  to  permit  the  insertion  of  voltage  calibra¬ 
tions  from  a  high  stability  source. 

Dynamics  measurements  in  the  re-entry  vehicles  have  been  of  three  types, 
rate  gyros,  differential  pressure  gauges  and  accelerometers.  Three  rate  gyros 
are  usually  installed  measuring  pitch,  yaw  and  roll  rate.  Power  supply  is  from 
a  rotary  invertor  and  an  inductive  bridge  type  of  pick-off  is  used  to  provide 
a  voltage  output  proportional  to  angular  velocity. 

The  acceleration  and  pressure  transducers  are  also  of  the  inductive  bridge 
type  but  power  requirements  are  modest  so  that  a  typical  installation  of  7 
transducers  can  be  small  ar.d  compact.  Such  a  system  would  consist  of  two 
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pressure  transducers  for  differential  pressure  measurements  in  the  pitch  and 
yaw  planes;  two  longitudinal  accelerometers,  one  for  acceleration  from  boosts 
prior  to  re-entry  and  one  for  deceleration  during  re-entry;  and  three  lateral 
accelerometers,  one  of  which  is  of  greater  sensitivity  to  monitor  motion 
during  the  early  part  of  re-entry.  All  these  transducers  are  capable  of  a 
full  scale  accuracy  under  the  most  severe  environment. 

In  most  Black  Knight  re-entry  vehicles,  where  many  other  measurements 
are  made,  there  is  insufficient  room  for  a  complete  dynamics  installation 
and  only  accelerometers  are  installed.  Such  a  system  is  illustrated  in 
Figs. 2  and  3«  Fig. 2  shows  the  sub-miniature  tape  recorder  and  its  armoured 
cassette.  Fig. 3  shows  an  instrumentation  assembly  with  the  tape  recorder 
visible  on  the  left  of  the  chassis.  As  we  have  seen  this  can  be  a  most 
useful  system  although  there  has  often  been  some  difficulty  in  placing  the 
accelerometers  in  the  optimum  positions. 

After  flight  the  tape  is  recovered  and  replayed  on  to  paper  records 
using  an  oscillograph  or  galva  recorder.  Timing  pulses  are  counted  down  to 
give  pulses  at  one  millisecond  intervals.  Synchronisation  with  ground 
instrumentation  is  achieved  by  simultaneous  recording  of  some  event  such  as 
a  boost  ignition.  The  paper  records  are  read  by  hand  using  a  semi-automatic 
method  at  intervals  of  0*01  seconds  or  less.  It  is  then  processed  using  a 
'Pegasus’  computer  to  produce  punched  paper  tape  of  the  parameters  as  a 
function  of  time.  A  print  of  this  may  be  plotted  for  the  purposes  of 
graphical  analysis  or  analysed  on  the  larger  'Mercury'  computer  by  the 
automatic  method. 

k  METHODS  OF  ANALYSIS 

The  automatic  method  of  analysis  is  essentially  a  curve  fitting 
technique.  Its  one  drawback  is  the  need  to  make  a  rough  guess  at  the 
parameters  of  the  curve  before  the  automatic  processes  can  proceed.  This 
is  because  the  dynamic  modol  of  re-entry  behaviour  is  a  highly  non-linear 
function  of  the  parameters  to  be  determined.  V/hereas  the  familiar  polynomial 
curve  is  a  linear  function  of  its  coefficients  and  may  bo  fitted  by  a  simpler 
and  more  direct  approach.  To  obtain  approximate  values  of  the  parameters  of 
the  dynamic  model  it  is  necossary  to  do  the  analysis  of  the  first  half  second 
or  so  of  ro-ontry  by  hand.  It  will  be  appropriate,  therefore,  to  briefly 
consider  such  an  essential  preliminary  to  the  automatic  technique.  At  the 
same  time  some  of  the  limitations  of  graphical  analysis  will  bo  demonstrated. 
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4.1  Graphical  methods 

The  production  of  a  trajectory  by  graphical  means  is  fairly  obvious. 

One  proceeds  from  drawing  a  curve  through  a  plot  of  height  against  time  data 
to  obtain  estimates  of  velocity  from  the  slope  and  thus  estimates  of  dynamic 
pressure.  A  summary  of  the  procedure  is  given  in  Appendix  A. 

Graphical  analysis  of  accelerometer  or  rate  gyro  data  to  determine  the 

motion  of  the  body  is  also  fairly  straightforward  and  has  been  described  in 

2  3 

many  papers  including  those  by  Nelson  *  .  An  illustration  is  given  in  Fig.4 
of  the  basic  idea  of  the  method.  It  shows  how  arm  lengths  R  and  S  and  fre¬ 
quencies  may  be  extracted  from  suitable  plots  of  the  data.  Applying  the  theory 
already  given  and  using  the  results  of  trajectory  analysis,  it  is  possible  to 

estimate  the  aerodynamic  derivative  C„  . 

Ma 

This  is  as  far  as  is  necessary  to  go  if  numbers  for  the  priming  of  the 
automatic  method  are  required.  All  t'at  is  needed  are  estimates  of  C^, 
p,  R  ,  S  and  the  phase  angles  y_  and  vc  at  the  initial  time  t  .  If  the 

0  0  ft  b  0 

o  o 

whole  analysis  is  to  be  conducted  by  hand  then  the  process  has  to  be  repeated 
for  convenient  lengths  of  the  whole  record  of  re-entry,  perhaps  lasting  a 
tota^.  of  10  seconds  or  more.  It  is  then  possible  to  obtain  plots  of  R  and  S 
against  time  and  to  make  estimates  of  and  by  studying  the  decay  of 
R  and  S  with  time.  Details  of  the  analysis  of  dynamic  motion  by  hand  are  also 
given  in  Appendix  A. 

There  would,  of  course,  be  no  need  for  an  elaborate  automatic  technique 
if  the  graphical  method  were  as  easy  as  it  sounds.  In  practice  data  does  not 
plot  with  tbo  beautiful  precision  of  the  smoothed  data  of  Fig.4  and  corrections 
for  the  effects  of  accelerometer  displacements,  for  instance,  are  difficult  and 
tedious  to  make.  Consequently  errors  tend  to  be  large  and  unknown  so  that  it 
is  often  impossible  to  be  sure  whether  irregularities  in  re-entry  behaviour  are 
real  or  apparent.  On  the  other  hand,  an  automatic  method  is  not  only  able  to 
remove  the  tedium  of  analysis  but  make  more  elaborate  corrections  for  known 
errors  and  treat  random  errors  in  a  proper  statistical  manner. 

4.2  The  automatic  method 

The  general  theory  of  the  automatic  method  is  given  in  Appendix  B.  It 
will  be  helpful  here  to  illustrate  Ihe  main  results  by  taking  the  simple  one 
dimensional  problem  of  fitting  a  trajectory  to  observation  of  height  at 
times  t^.  We  will  suppose  that  heights  can  be  calculated  from  the  function 


RESTRICTED 


18 


RESTRICTED 


*ic  s  *(V  V  V 


(26) 


and  it  is  desired  to  determine  x^  and  x^,  which  could  be  initial  velocity 


and  height,  so  as  to  minimise  the  function 

U  =  2  wi(*i  -  *iQ)2 

which  is  the  sum  of  the  squares  of  the  weighted  residuals. 


(27) 


Because  %  is  a  non-linear  function  of  x^  and  x we  employ  the  method 
of  differential  corrections  which  starts  with  approximate  values  x^  for 
and  obtains  x^  by  an  iterative  procedure.  That  is,  we  start  with 

(28) 


ic 


and 


V 

=  xk  +  6xk  (k  =  1 ,  2)  . 


(29) 


The  technique  is  to  expand  equation  (28)  by  Taylor's  series  in  terms  of  the 
partial  derivatives 

a*(v  V 


ik 


It  is  shown  in  Appendix  B  that  Sx^.  is  given  by  solving  the  equation: 


¥  E  a  CWD 


(30) 


whore 


/Sxl\  •  C  s  /fil,,,fnl\  » 
\  x2  /  \i’l2.“fn2/ 


*io  “  *i 


\ 


and  Y  =  CW  C* 


,%  -  *  y 

^  no  n  ' 

where  asterisk  *  denotes  the  transpose.  In  long-hand  (30)  becomes: 

6x1  2  w1fi1fi1  +  6x2  2  w1fi1f12  ■  s  wifii^io"xi^ 

6x1  1  Vi1fi2  +  5x2  2  Vi2fi2  "  2  Vitale  "*i)  (32) 
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these  are  known  as  the  Normal  Equations  and  are  solved  by  a  pivotal  condensa¬ 
tion  scheme  to  obtain  6x^.  The  non-linear  nature  of  the  problem  makes  it 
necessary  to  repeat  the  process  using  the  new  values  of  until  6x^  is 
negligible. 

The  accuracy  of  the  observations  ia  estimated  from: 


2 


cr 


2  E*  CWD 

cr - 

o  n-  p 


where  cr*  is  the  variance  of  the  uncorreoted  weighted  observations. 

At  the  95/'a  probability  level,  the  probable  error  in  the  parameters  x^. 
is  given  by 


2  -1 

AX  =  2  cr  «k 

-1  -1 
where  i3  the  kth  diagonal  element  of  the  matrix  V  . 

The  partial  derivatives  are  computed  by  a  simple  perturbation  technique. 
In  thi3  case  three  trajectories  would  be  calculated  instead  of  only  one,  tho 
first  for  ,  x^,  the  second  for  (x^  +  ^)»  x2>  and  the  third  for  x^ ,  (x2  +  £2). 
Hie  partial  derivatives  are  then  given  by: 

x(x1 ,  x2,  t±)  -  x(x1  +  ,  x2,  t^ 

fn  =  ej 

and  similarly  for  f^. 

Hie  purpose  of  the  weight  is  to  give  the  most  weight  to  the  most 
accurate  observations.  They  are  therefore  defined  as  numbers  inversely 
proportional  to  the  square  of  tho  probable  errors  of  tho  observations,  i.e.  the 
reciprocal  of  the  variances  of  the  individual  observations 


wi  = 


In  practice  it  is  found  convenient  to  read  in  (l/6^)  to  tho  computer,  rather 
than  tho  actual  weight. 

When  data  is  weighted  in  this  way,  the  estimate  of  cr,  the  weighted 
standard  deviation  of  tho  residuals,  should  bo  little  more  than  unity.  Values 
of  2  or  3  for  cr  would  indicate  that  6^  aro  unduly  optimistic  or  that  there  are 
defects  in  the  dynamic  model.  For  most  data  6^  aro  not  available  in  which  case 
all  the  observations  are  given  unit  weight. 
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This  then,  very  briefly,  is  tho  theoiy  of  the  automatic  method.  Although 
a  simple  case  has  been  taken  the  theory  is  completely  general  and  has  been 
programmed  for  a  computer  for  any  number  of  parameters,  dimensions  and 
observations.  Into  this  general  programme  the  subroutine,  for  tho  calcula¬ 
tion  of  the  dynamic  model  appropriate  to  the  observations,  is  inserted. 

A  block  schematic  diagram  of  the  complete  programme  is  shown  in  Pig. 5  and 
details  are  given  in  Appendix  C.  Tho  main  features  are  as  follows: 

(a)  Data  from  the  complete  re-entry  period  is  split  into  blocks  of  a 
length  depending  on  how  long  parameters  may  be  considered  constant.  It  is 
only  necessary  to  estimate  initial  conditions  for  the  first  block,  the  pro¬ 
gramme  predicts  for  the  remainder. 

(b)  A  facility  is  incorporated  for  rejecting  data  points  which  are 
obviously  wrong  as  a  rosult  of  misreading;  this  is  on  the  basis  of  tho 

residual  being  greater  than  4  times  the  rms  of  the  whole.  Only  data  accepted 

2 

is  counted  and  its  variance  o'  calculated.  The  programme  fails  if  data 
accepted  is  less  than  p,  tho  number  of  parameters.  This  usually  arises  if 
the  solution  diverges. 

(c)  The  fitting  procedure  is  applied  several  times  to  each  block 
until  tho  values  of  6x^  are  less  than  a  preset  accuracy  level.  The  programme 
then  goes  on  to  the  next  block  after  printing  out  extra  information  such  as 
the  fitted  trajectory  or  motion  and  the  rosiduals. 

5  AUTOMATIC  ANALYSIS  RESULTS 

5.1  Trajectory  analysis 

Results  to  be  quoted  wore  obtainod  by  a  relatively  simplo  one  dimensional 
fit  to  height  data,  using  the  model  discussed  in  Section  2.  Those  particular 
programmes  form  only  an  appendix  to  a  complete  system  for  providing  smoothed 
trajectories  from  launch  to  impaot  in  three  dimensions  for  the  whole  vehicle. 
However,  re-entry  test  vehicles  usually  follow  noar  vertical  trajectories, 
so  that  one  dimensional  fits  to  height  data,  with  an  allowance  for  the 
departure  of  the  trajectory  from  tho  vertical,  have  been  found  adequato  for 
re-entry  dynamics  studies. 

Two  programmes  aro  in  conmon  uso  -  TRAJECTORY  ANALYSIS  1623/lA  and 
1623/lB  -  and  fit  a  trajectory  to  height  data  derived  from  ballistic  cameras 
or  radar.  Both  programmes  use  atmosphoric  pressures  and  temperatures  measured 
at  the  time  of  flight.  Programme  1A  uses  a  drag  table  of  C^  against  Maoh 
number  based  on  the  best  pro-flight  estimato.  The  first  parameter  is  allocated 
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to  AC^,  the  mean  zero  error  in  the  drag  table.  Second,  and  ttiird  parametors  are 

initial  height  h  and  velocity  v  .  No  attempt  is  made  to  allow  for  induced 

drag  from  the  dynamic  motion  of  the  vehicle.  Programme  IB  uses  measured  values 

of  acceleration  A  and  A  in  place  of  the  drag  table,  otherwise  it  is  similar 
x  y 

to  1A.  The  first  parameter  is  now  AAx>  the  zero  error  in  longitudinal 
acceleration.  The  dynamic  model  is  given  by  equation  (23)  et  seq  and  makes 
full  allowance  for  induced  drag  and  accelerometer  displacement  effects.  The 
programme  calculates  the  drag  coefficient  C^.  Operating  details  of  these 
programmes  will  be  found  in  Appendix  C.  The  form  of  a  typical  set  of  data 
for  input  to  the  computer  being  reproduced  in  Table  1. 

Ihe  result  from  Programme  IB  is  reproduced  in  Table  2.  It  will  be  noted 
that  the  programme  first  prints  the  initial  parameters  as  read  in,  the  seoond 
column  being  the  desired  perturbations  Subsequently  the  new  values  of  the 

parameters  are  printed  after  each  iteration  but  the  second  column  now  lists  the 
best  estimates  of  the  standard  deviations  AX  of  the  parameters.  Also  before 
each  new  parameter  list  the  best  estimate  of  Ihe  standard  deviation  of  the 
weighted  residuals  is  printed,  followed  by  the  number  of  degrees  of  freedom 
(n-p)  where  n  is  number  of  observations  accepted. 

It  will  be  noted  that  one  iteration  sufficed  to  obtain  the  best  values 
of  the  parameters,  the  seoond  made  no  further  improvement.  In  the  printed 
trajectory,  the  unweighted  residuals  are  printed  on  the  right  hand  side.  The 
initial  fall  in  with  decrease  in  H  is  of  particular  interest  *uvd  has  been 
ascribed  to  a  viscous  effect.  Only  part  of  the  fitted  trajectory  is  reproduced 
for  reasons  of  brevity. 

5.2  Dynamic  analysis 

Four  programmes  havo  been  developed  at  present  for  the  successful  analysis 
of  flight  data,  one  each  for  differential  pressures  and  rate  gyro  analysis  and 
two  for  accelerometers,  as  follows: 

Dynamic  analysis 

1623/1  “  For  differential  pressures  - 
9  parameters 

1623/2  -  For  accelerometers  -  9  parameters 
1623/3  ”  For  acoelorometers  -  10  parameters 
1623/4  “  For  rate  gyros  -  9  parameters 

The  extra  parameter  in  the  case  of  No. 3  is  C^,  the  other  9  aro  C^,  C^,  P, 

R  ,  S  »  Ye  and  two  zero  errors  E  ,  B  .  At  present  programmes  1  to  3 

o  0  k  o  x  y 
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cannot  distinguish  between  trim  and  zero  errors  but  later  versions)  it  is 
hoped,  will  do  so. 

All  programmes  use  the  same  subsidiary  data  as  Trajectory  Analysis  1A 
exoept  that  the  drag  table  is  derived  from  IB  results,  as  are  the  other 
trajectory  parameters.  The  trajeotory  is  needed  to  give  values  for  density 
and  dynamic  pressure  and  is  obtained  by  integration  of  equations  (22) ,  (24) 
and  (25).  Allowance  is  made  for  induced  drag  so  that  the  oaloulated  trajeo¬ 
tory  is  identical  to  that  deduced  by  trajectory  analysis.  It  should  be  noted 
that  only  and  Cjja  are  read  into  the  machine,  is  calculated  from  the 
difference  (c^  *»  -  C^).  Thus,  exoept  in  the  case  of  Programme  No. 3, 

is  fixed  and  damping  can  only  be  altered  by  varying  C^. 

All  programmes  may  be  used  to  fit  either  a  single  instrument  or  an 
orthogonal  pair .  It  is  thus  possible  to  do  a  oomplete  analysis  from  one 
instrument  only.  Howaver,  Programme  3  is  intended  to  analyse  two  accelero¬ 
meters  at  different  axial  stations  for  the  determination  of  C^,  rather  than 
rely  on  differential  damping  of  R  and  S  to  give  and  thus  C^.  Operating 
details  aro  given  in  detail  in  Appendix  C. 

Up  to  date,  the  automatic  method  has  been  used  to  analyse  short  lengths 
of  reoord  or  blocks  of  data.  The  complete  record  is  coverod  by  a  series  of 
blocks.  Some  examples  of  the  analysis  of  a  single  blook  are  shown  in 
Tables  3  and  4  and  Figs. 6-8.  Table  3  is  the  data  tape  and  Table  4  the  result 
of  fitting  data  from  a  pair  of  aooelerometers  with  Programme  No. 2.  The  result 
is  plotted  in  Fig.6.  Actual  measurements  are  plotted  for  comparison  with  the 
smoothed  curve  of  aooele ration.  The  motion  plotted  in  Fig.4  to  illustrate 
graphical  analysis  was  obtained  from  the  same  accelerometer  analysis  of  Fig.6 
and  is  thus  the  required  end  produot  of  the  analysis,  it  is,  of  course, 
dependent  on  the  assumed  value  for  C^.  It  will  be  noted  that  the  pattern  of 
Fig.4  and  Fig.6  is  the  result  of  R  and  S  being  almost  equsl,  that  of  Fig. 7 
and  Fig. 8  is  for  a  case  when  2R  a  S.  Programme  2  was  again  used  to  analyse 
a  pair  of  aooelerometers  to  give  Fig. 7.  Pro grannie  4  was  used  to  analyse  a 
rate  gyro  record  from  the  same  re-entry  but  at  on  earlier  time  to  give  Fig. 8, 
Some  of  the  parameters  obtained  from  this  analysis  are  given  along  side  those 
from  an  accelerometer  analysis  for  the  same  part  of  the  re-entry.  It  will  be 
noted  that  values  of  the  parameters  agree  very  well  but  the  rate  gyro  is 
claimed  to  be  much  more  acourato.  This  is  probably  beoause  a  gyro  is  not  at 
the  mercy  of  so  many  side  effects  as  an  aooelerumeter. 
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Tha  poor  accuracy  of  Cu  in  the  figures  quoted,  in  Pig. 8,  especially  in 

Mq 


the  case  of  the  accelerometer,  is  mainly  because  the  record  is  for  an  altitude 
of  160  000  ft  when  the  effect  of  Cy^  is  negligible.  Accuracy  is  greatly 
improved  at  lower  altitudes  and  is  also  helped  by  analysing  a  longer  period  of 
record.  Even  so,  Cy^  is  a  difficult  parameter  to  measure  accurately  and  it  is 
doubtful  if  an  accuracy  of  bettor  than  ±0*5  can  bo  achieved. 


It  has  been  seen  that  trajectory  analysis  converges  very  rapidly. 

Dynamic  analysis  is  much  slower,  in  spite  of  good  guesses  at  the  initial 
parameters,  because  of  the  highly  non-linear  form  of  the  dynamic  model.  This 
is  illustrated  by  the  plot  of  parameter  values  against  iteration  number  in 
Fig. 9  for  an  accelerometer  analysis.  The  standard  deviation  of  the  residuals 
and  the  number  of  degrees  of  freedom  are  also  plotted,  and  it  will  be  seen 
that  initially  half  the  observations  were  not  used.  This  probably  accounts  for 
some  of  the  erratic  variations  in  parameter  values  as  more  observations  were 
accepted.  It  has  been  found  in  practice  that  up  to  six  iterations  are 
nocessary  to  fit  the  motion  at  the  beginning  of  re-entry  when  the  initial 
estimates  are  the  result  of  rough  graphical  analysis.  At  present  a  complete 
re-entry  record  is  split  into  10  blocks  of  about  AO-50  readings  and  covering 
20000  ft  of  altitude.  Having  made  an  initial  guess  at  the  first  block  the 
remaining  nine  are  handled  automatically  in  about  one  hour  of  computing.  The 
time  is  the  same  whether  one  instrument  or  an  orthogonal  pair  are  boing 
analysed. 


6  FUTURE  DEVELOPMENTS 


This  paper  has  been  primarily  concerned  with  the  basic  approach  to  the 
analysis  of  re-entry  data  by  automatic  moans.  At  present  this  soheme  uses  the 
simplest  ctynamic  model  which  can  reasonably  bo  assumed  -  the  caso  of  linear 
aerodynamics  and  a  symmetric  missile.  It  has  been  found,  however,  that  it 
works  surprisingly  well  and  is  adequate  for  a  primary  analysis  of  re-entry 
data. 


The  analysis  of  several  firings  on  this  basis  has  shown  that  thoro  are 
at  least  two  defects  in  the  model  which  must  be  corrected  boforo  analysis  of 
those  firings  can  bo  satisfactorily  completed.  The  most  common  defect  is  the 
non-linear  behaviour  of  Cyft.  It  is  this  which  makos  it  nocessary  to  analyse 
over  unduly  short  time  intervals.  To  analyse  over  longor  poriods  assuming  a 
constant  Cya  is  impossible  booause  governs  the  aerodynamic  frequenoy  «  and 
small  errors  soon  build  up  into  unacceptably  large  phase  errors.  As  a  .a suit 
of  taking  data  in  short  groups  it  is  impossible  to  obtain  a  smoothed  ourve  of 
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incidence  through  re-entry  and  accurate  estimates  of  C 
essential  to  permit  to  vary  in  some  way. 


wq 


and  C.  . 
La 


It  is  thus 


Results  have  shown  that  incidence  has  the  dominant  influence  on  and 
that  a  function  of  the  form 


C 


Ua 


-2 

c  +  o4  a 
0  1 


(36) 


is  a  good  approximation  to  the  truth.  It  is  proposed  to  extend  the  dynamic 
model  in  this  way,  at  the  same  time  putting  a  restriction  on  spin  rate  P. 

This  is  necessary  because  non-linearities  in  also  manifest  themselves  as 
apparent  spin.  Spin  will  bo  set  at  the  best  known  free  spaoe  value  if 
nothing  better  is  available.  This  modification  should  enable  the  length  of 
record  that  may  be  analysed  to  be  greatly  extended. 

It  will  bo  noted  that  Cu  will  be  made  a  function  of  the  mean  incidence 
—  2  2 

a  or  VR  +  S  ,  however,  it  it  is  really  a  function  of  a  -ip  and  not  only 
affects  pitching  frequency  but  the  waveform  of  the  oscillations.  This  ought 
to  be  taken  into  aocount.  Theoretical  work  on  this  problem  by  Rasmussen^  and 
others  has  been  successfully  applied  to  ballistic  range  experiments.  Nelson^ 
has  extended  the  non-linear  theory  to  a  varying  environment  on  an  intuitive 
basis.  All  these  workers  use  graphical  methods  of  analysis  and  it  is  not 
clear  how  to  apply  these  methods  in  an  automatic  scheme.  Moreover  it  is 
doubtful  whether  re-entry  data  is  of  sufficient  quality  to  Justify  a  detailed 
non-linear  analysis  soheme. 


A  further  limitation  of  the  dynamio  model  is  the  assumption  of  aero¬ 
dynamic  symmetry.  On  at  least  one  round  this  defect  has  proved  serious  and 
has  made  it  impossible  to  fit  data  over  long  periods  up  to  resonance,  because 
of  the  shift  in  phase  and  amplitude  of  trim.  After  resonance  the  trim  is  in 
phase  with  the  S  arm  and  gives  no  further  trouble.  Equation  (4)  in  Section  2 
includes  the  effect  of  asymmetry  and  it  was  mentioned  that  the  solution  is 
the  sum  of  two  particular  integrals.  That  for  the  motion  about  trim  forms 
the  present  dynamic  model,  it  remains  to  add  the  particular  integral  for  the 
behaviour  of  trim.  An  analytical  solution  for  this  has  been  obtained  by 
G.S.  Green  at  the  R.A.E.  for  tho  case  of  no  damping  and  lift,  it  is  a  series 
arising  from  a  modified  fora  of  Bessel's  Function.  The  series  converges 
rapidly  well  aw^r  from  resonance  and  is  useful  for  determining  the  general 
effeot  of  asymmetry.  In  the  region  of  resonance  it  oan  go  to  hundreds  of 
terms  and  is  quite  unsuitable  for  our  purpose.  Some  other  approach  is 
desired. 
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It  has  been  decided  that  the  simplest  way  to  determine  the  particular 
integral  for  the  motion  of  trim  is  to  integrate  equation  (b)  numerically  with 
a  -  i^  zero  at  infinity.  The  motion  about  trim  would  still  use  the  analytical 
solution.  This  is  not  as  bad  as  it  sounds  because  ae-i£eis  a  scale  factor  and 
the  only  important  variables  are  and  perhaps  Cy^.  It  would  involve  adding 
12  differential  equations  to  the  mathematics  of  the  programme  and  should  oauso 
no  undue  difficulties. 

Prom  experience  with  the  automatic  scheme  using  the  linear  symmetric 
model  it  is  felt  that  the  above  two  modifications  will  enable  a  complete 
re-entry  to  be  analysed  in  a  single  run  on  the  computer.  ftiis  will  not 
necessarily  speed  the  analysis  up  but  will  enable  accurate  mean  values  for 
Cyo»  CT^  and  Cy^  to  be  obtained  as  well  as  a  continuous  history  of  incidence 
with  time.  It  is  hoped  to  try  it  out  in  the  first  four  months  of  1965* 

7  CONCLUSIONS 

The  automatic  method  of  analysis  has  proved  to  be  a  considerable  advance 
on  the  graphical  methods  formerly  used  for  the  analysis  of  re-entry  data.  Not 
only  has  the  primary  objective,  of  developing  a  quick  and  effioient  method  of 
treating  large  quantities  of  re-entry  data,  boon  achieved,  but  it  has  enabled 
a  more  sophisticated  instrumentation  scheme  using  accelorometers  to  be 
recommended.  This  was  an  unexpected  bonus  and  moans  that  a  very  simple 
arrangement  of  accelerometers  can  monitor,  not  only  the  motion,  but  all  the 
aerodynamic  characteristics  of  the  re-entry  vehiole.  This  principle  might  be 
applied  with  profit  in  ballistic  range  experiments. 

Computer  methods  also  enable  more  complex  mathematical  models  of  motion 
to  bo  used  than  is  normally  practicable  in  hand  analysis,  such  as  the  treatment 
of  non-linear  aerodynamics  and  asymmetric  effects.  In  addition  experimental 
errors  are  treated  in  a  consistent  manner  which  is  important  when  results  of 
several  firings  are  to  be  comparod.  Those  considerations  are  quite  general  and 
do  not  only  apply  to  re-ontry  work.  It  is  probable  that  such  methods  can  be 
applied  to  other  experiments  in  the  aero-space  field.  Re-entry  experiments 
cannot  be  the  only  case  where  simplor  instrumentation  can  be  gained  at  the 
expense  of  more  elegant  analysis  methods. 
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Appendix  A 

GRAPHICAL  METHODS  OF  ANALYSIS 


The  term  graphical  methods  includes  all  the  many  techniques  for  the 
analysis  of  flight  data  which  do  not  require  a  digital  computer;  although 
it  is  difficult  to  draw  a  definite  dividing  line.  The  ideal  automatic 
method  should  not  need  the  intervention  of  a  human  operator  but  all  methods 
of  the  curve  fitting  type  roquiro  a  guess  at  tho  initial  parameters  and  it 
is  here  that  graphical  methods  must  be  employed.  Graphical  mothods  therefore 
form  an  essential  preliminary  to  an  automatic  analysis  and  have  not  been 
included  for  historical  interest.  These  techniques  vary  greatly  with  the 
data  available  and  the  whim  of  the  analyst  and  are  extensively  described  in 
tho  literature.  Only  a  brief  outlino  will  be  given  in  this  paper. 

Trajectory  analysis 

The  stops  in  the  process  would  normally  include  the  following: 

(a)  Tho  synchronisation  of  ground  and  airborne  time  which  can  bo  tricky 
when  an  airborne  tape  recorder  is  in  use.  The  usual  method  is  to  record  a 
flash  both  on  tho  ground  and  on  board  at  some  convenient  time. 

(b)  The  determination  of  velocity  V  and  height  H  at  a  time  shortly 
bofore  re-entry  (say  250  000  ft  altitude).  This  will  be  provided  by  analysis 
of  vacuum  trajectory  data  such  as  from  radar  and  is  a  subject  in  itself. 

(c)  The  integration  of  the  longitudinal  acceleration  record,  starting 
with  the  initial  V  and  H,  to  obtain  velocity  and  height  as  functions  of  timo. 

2 

(d)  The  calculation  of  -g  pV  using  measured  pressures  and  temporaturos 
of  the  atmosphere  at  tho  time  of  firing. 

(e)  A  cross  chock  and  adjustment  (if  necessary)  of  the  trajectory  with 
spot  heights  from  Ballistic  Camera  Data. 

(f)  The  estimation  of  tho  drag  coefficient. 

Dynamio  analysis 

Tho  analysis  would  include  some  or  all  of  tho  following  steps  although 
not  necessarily  in  the  order  given.  Tho  components  of  the  record  from  a 
particular  instrument  are  referred  to  as  R'  and  S’ ,  and  could  be  actual 
incidence  or  angular  volocity  or  some  othor  variable. 


RESTRICTED 


38 


Appendix  A 


RESTRICTED 


27 


(a)  For  tho  oarly  part  of  ro-ontry  whon  pitching  rate  u>  is  low  and  of  tho 
same  order  as  spin  rate  P,  oxtract  tho  froquonoy  and  amplitude  of  R*  and  S'  by 
study  of  a  record  from  a  single  instrument.  Fig. A  demonstrates  in  diagrammatic 
form  how  this  can  be  done.  By  application  of  equation  (6)  «  and  P  may  be 
extracted  and  plottod  against  time,  along  with  R'  and  S'.  Ibis  typo  of 
approach  for  just  tho  first  few  re-entry  oscillations  is  all  that  is  required 
to  prime  the  automatic  system. 

(b)  As  an  alternative  to  tho  above,  especially  whon  w  is  large,  plot  tho 

pitch  readings  against  yaw  readings  to  obtain  a  graph  similar-  to  Fig. 4.  Hve 

precise  form  of  this  will  depend  on  tho  relative  sizes  of  w  and  P,  and  R'  and 
2 

S'  .  Measure  the  time  and  angle  between  adjacent  peaks  as  shown,  perhaps 
taking  several  peaks  to  obtain  a  mean  value  for  a  chosen  time.  Note  that 
angles  must  bo  measured  from  tho  trim  centre. 

Then 

<u  s  V(t2"'tp  and  P(l  -  A./2B)  =  (it- A0)/(t2-  t^) 

where  t2  and  t^  corresponds  to  time  between  peaks. 

Measure  the  amplitude  of  maxima  and  minima  from  the  trim  centre  as  shown. 
The  maximum  equals  (R’ +  S')  and  the  minimum  iB  (R'-  S')  so  that  R'  and  S'  may 
be  obtained  as  before. 

(o)  A  plot  of  R'  and  S'  against  time  is  obtained  for  each  typo  of 
instrument  or  different  pair?  of  accelerometers  by  repeated  application  of  (a) 
or  (b)  abovo  at  a  number  of  different  timos.  From  u>,  P  and  £ pV  from  tho 
trajectory,  R'  and  S'  can  be  procossed  to  obtain  the  incidence  components  R 
and  S,  the  pitching  moment  derivative  and  perhaps  the  normal  force 
derivative  C^.  Ibis  is  achiovod  by  application  of  tho  appropriate  equations 
in  Sections  2  and  3. 

(d)  Tho  final  stage  in  tho  analysis  is  tho  stuty  of  the  variation  of  R 

and  S  with  time  to  derive  the  stability  derivatives  CT  and  C  .  First  it  is 

ua  mq 

necessary  to  remove  tho  effect  of  the  term  -  often  the  dominant  torm  for 

much  of  tho  re-entry.  Hence  compute  and  plot  against  time: 

R  *  /ffi  R  and  S  s  Vo5  S  • 

Tho  slopes  of  tho  resulting  corrected  curvos  diydt  Bnd  d2/dt  are  then  measured. 
Then  it  follows  from  equation  5 
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a  +  -*  2B 


-  2B 
u 


*  at  *P5 


Since  a  and  b  are  functions  of  and  Cj^  it  is  theoretically  possible  to 
solve  these  equations  to  obtain  the  derivatives.  Accuracy  i3  usually  too 
poor  to  make  this  possible  but  can  be  obtained  separately  from  knowledge 


of  C„  and  C,  for 
Na  A 


-  C_.  -  C. 
Na  A 


Hence  a  mean  value  for  C„  may  bo  determined. 

Nq 
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GENERAL  THEORY  OF  AUTOMATIC  METHOD  OF  ANALYSIS 

The  differential  correction  procedure) 

In  general,  a  dynar.ic  observation  will  be  a  vector  such  as  position, 
velocity,  or  acceleration,  and  is  best  defined  by  the  components  in  some 
cartesian  axes  system.  In  practice  only  one  or  two  particular  components 
may  bo  measured,  but  the  following  theory  assumes  three  dimensions  although 
it  is  equally  valid  for  any  number. 

Let  a^,  |3^,  b®  observations  of  the  vector  at  a  timo  t^,  and  let  n 
such  observations  be  made  at  different  times.  If  there  are  P  undetermined 
parameters  (k  a  1 ,  2,  ...  p)  in  the  dynamic  model,  and  if  there  are  no 
errors  in  the  observations  or  the  model  then  ti.o  observations  would  bo  related 
to  the  parameters  by  equations  of  the  form 

-  a(X^ ,  X2,  ...  Xp,  t^) 

•  P(V  x2,  •••  Xp.  V 

-  Y^,  Xg,  ...  Xp,  t^) 

the  functions  a,  0,  y  depending  on  the  dynamic  model.  VJhen  Jn  >  P,  the 
parameters  may  be  obtained  by  solving  tho  equations  redundantly,  using  a  least 
squares  procedure  to  determine  tho  best  values. 

Let  a^Q,  y^c  b®  computed  values  close  to  the  observed  values  based 
on  estimates  of  tho  parameters  x^,  such  that 

aic  e  a(v  v  —  v  V 

Pic  s  P(V  X2*  •••  v  V 

Yi0  =  Y(^,  x2,  ...  Xp,  t±) 

and  let  tho  residuals  bo  of  tho  form: 


R  .  =  a.  -  a,  etc  . 

ai  i  ic 


Then  the  procedure  adopted  is  to  determine  tho  values  of  x^  so  as  to  minimise 
the  quantity 
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0  ■  Y  V  *«i  *  l,  V  EPi * 


inO 


n_  n_ 

Y  v  &  4  E  v *• 

i=0  i=0 


2 

LYi 


where  w  wa.  and  w  are  the  weights  of  the  observations.  Suppose  are 
ai  pi  yi 
approximations  to  x^  so  that 

*k  "  *k  *  6\ 

•1.  •  a(xi  ■  V  •  ete-  <B,) 

Linear  equations  may  be  obtained  by  expanding  equation  (Bl)  as  a 
Taylor  series ,  thus 


“io  *  “io  *  £  fkl  ***  4  0C(6’‘k)2  ’  *l> 
k=1 
P 

K'lc  ’  he  *  Y  8ki  4  0(<8*k>2  ’  ^ 


k=1 

P 


Yic  =  Yio  +  E  hki  ***  +  °((6Xk)2  ’  ^ 


where 


da 


fki  ° 


k=1 


t,)  dp(Xj:,  t.t) 

- ’  gki  = 


exk 


*v(y  V 

\i  "  dxk 


2 

The  correction  procedure  is  based  on  the  assumption  that  (&x^)  ^-s 

negligible.  In  practice  this  is  not  generally  true,  and  the  minimisation  of 
U  is  obtained  by  repeated  application  of  the  procedure. 


Lot 


Ki  ’  *i  •  “lo  ’  oto 


so  that 
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and 


* 

Rdi  =  Rai  "  fki  6xk 
k=1 

Rpi  =  V  “  5_,  8ki  6xk 

k=1 
P 


a  ,  =  n.  . 

fi  Yi 


'Y,  hkitait 


k=1 


(B2) 


^  (Wai  Rdi2  +  V  Rpi2  +  Wyi  Rdi2)  * 
i=1 


(B3) 


It  is  required  to  determine  Sx^  so  ns  to  minimise  U’,  that  is  to  find 
values  of  for  which  dU'/d^x^)  is  zero.  For  this  purpose  it  is  helpful  to 
pose  the  problem  in  matrix  form. 


Let 


c  « 

/f11 ' 

f12  »  *  ’ 

g11 

*  *  *  » 

h11 

»  h12 

f21  ’ 

f22  »  *  » 

g21 

*  *'22  *  *  ’ 

h21 

*  h22 

•  t 

l 

*  1  *  > 

• 

i  •  >  •  » 

• 

»  • 

'  fk2  ’ 

fk2  ’  *  » 

Sk1 

*  gk2  »  ’  » 

*  \2 

E  SS  /ox 


1 

6V 


and  W  a  diagonal  matrix  with  the  weights  along  the  diagonal  in  the  order 
(wai,  wa2»*,,wpi*  wp2’’*,WYl'  Wy2^  °^i0r  ^QrmB  bsins  *oro.  So  that 
equations  (B2)  and  (B3)  may  be  written: 

D'  =  D  -  C*  E 
U'  =  D'*  W  D' 

whore  the  asterisk  *  denotes  the  transpose. 


(B4) 

(B5) 
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From  equations  (B4)  and  (Bj) 


D'*f  D1  =  (D*  -  E*  C)  W(D  -  C*  E)  (B6) 

therefore 


3U»  _  a(D«»WDl) 

3(6^)  =  d(6xk) 


=  -  2  CWD  +  2  CWC*  E  =  0 


that  is 


CWC*  E  =  CWD 


(B7) 


which  is  the  matrix  form  of  what  are  commonly  called  the  normal  equations. 
Defining 


and 


¥  =  C  V.'  C* 

¥E  a  CWD 
E  =  CWD 


which  is  the  dosired  solution  and  gives  the  increment  to  be  added  to  the 
parameters.  The  procedure  is  repeated  using  the  now  parameters  as  starting 
values,  until  E  is  negligible  when  (Sx^  must  be  negligible,  and  the  best 
estimate  of  the  parameters  obtained. 

Estimation  of  the  accuraoy  of  the  observations 

To  estimate  tho  accuracy  of  the  observations  it  is  necessary  to 
determine  their  covariance  matrix.  Let  us  suppose  that  D  is  redefined  as  a 
voctor  of  residuals  based  on  the  true  values  of  the  parameters.  That  is  a 

,  \)). 

Define  D  as  tho  moan  value  of  D.  Then  the  weighted  covariance  matrix  of  D  is 

Cov(d)  =  Expectation  (D-D)  W(D*-D*)  .  (B8) 

To  determine  Cov(D)  cortain  assumptions  have  to  be  made:  that  the 
dynamical  model  exactly  fits  the  true  observations  and  that  tho  errors  are 
random,  unoorrolatod  and  normally  distributed.  In  which  oa3e  the  covariance 
matrix  reduces  to: 

Cov(d)  a  Expectation  (D*WD)  b  o ^  (B9) 

2 

a  best  estimate  of  a  is  obtained  from  tho  minimum  value  of  U*,  given  by 
substituting  (B7)  into  (B6). 


typical  oloment  in  D  is  [a^  -  aCX^,  t^)]  rather  than  [a^  -  aCx^ 
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Thus 

D'*WD*  »  D*WD  -  D*W  C*  E  -  E*  CWD  +  E*  CWC*  E 
»  D*WD  -  D*WC*  E  -  E*  ¥  B  +  B*  ¥  B 


Best  estimate  of 


so 


or 


D*W  D  -  E*  CWD 


cr  = 


U' 


3n-P 


(D^ffD1) 

3n-P 


D»W  D  E*  CWD 
3n  -  P  "  3n  -  P 


2  B*  CWD 
°o  "  3n-P 


(BIO) 


whore  a ■  is  the  variance  of  the  unoorreoted  weighted  observations. 
Estimation  of  the  aocuraoy  of  the  dvnamio  parameters 


An  estimation  of  the  aocuraoy  of  the  parameters  determined  by  the 
correction  procedure  is  most  important,  otherwise  there  will  be  a  danger  of 
drawing  unjustified  conclusions  from  the  results. 

Let  E  be  defined  as  a  vector  of  deviations  from  the  true  values  of  the 
parameters  X^.  ftiat  is  a  typical  element  of  fi  is 

covariance  matrix  of  E  is  given  by  taking  the  expectation  of  the  product 
(E-E)  (E*  -1*) . 


(X^-  x^) .  Then  the 


Now 

so  that 


E-E  =  ¥_1  CW(D-D) 


Exptn  (E-B)(E*-E*)  =  Exptn  ¥_i  C  W(D-5)(D-D)*  W*  C*  ¥~1# 
from  equation  (B9), 

Exptn  W(d-5)(D -5)*  »  cr2  . 

Hence 

Cov(B)  s  a2  ¥_1  CW*  C*  ¥-1* 

noting  that  ¥  and  W  are  symmetric  matrioos  so  that  they  are  unchanged  by 
transposition,  wo  have 
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Cov(E)  b 


(B11) 


Provided  the  parameters  are  unoorrolated  all  the  terms  in  Cov(K)  will  vanish 
exoept  for  the  diagonals.  These  diagonal  terms,  therefore  give  the  variance 
of  the  individual  parameters,  since  the  variances  of  6X,  must  be  the  variance 


of  Xk. 

Let  the  variance  of  X^  be 
f-1,  then 


sk  “4 


be  the  kth  diagonal  eloment  of 


(B12) 


The  best  estimate  of  the  standard  deviation  of  will  be 
correction  may  be  necessary  for  a  small  number  of  observations  (usually 
(n-P)  <  20).  It  is  usual  to  quote  accuracy  at  the  95$>  probability  level, 
in  which  case  the  probable  error  in  X^  is  given  by 

AX  «  2  Sk  =  2  o2  y”1  .  (B13) 

2 

It  should  be  emphasised  that  this  estimate  is  only  valid  when  (6x^.)  is 
negligible,  a  situation  which  is  usually  reached  after  several  repetitions 
of  the  correction  procedure.  The  non  diagonal  terms  in  i  ^  must  also  be 
negligible. 


S^,  although  a 
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OPERATING.  DETAILS  OF  THE  AUTOMATIC  ANALYSIS  PROGRAMMES 

The  automatic  scheme  described  in  Appendix  B  is  completely  general  and 
could  be  applied  to  all  manner  of  different  problems.  It  has  boon  developed 
to  meet  the  particular  needs  of  re-entry  data  analysis  and  advantage  has 
therefore  been  taken  of  features  peculiar  to  these  problems.  As  far  as 
possible  a  standard  procedure  has  boen  adopted  for  both  trajectory  and 
dynamio  analysis  programmes,  because  it  considerably  eases  the  task  of  writing 
new  programmes  and  modifying  others.  The  programmes  have  been  written  for 
exclusive  use  on  a  Mercury  Computer.  For  this  reason  much  of  tho  programme  is 
in  machine  code  mixed  with  Autocode,  with  the  aim  of  achieving  tho  maximum 
computing  speed.  This  would  limit  the  use  of  the  programmes  on  other  computers. 
Details  are  given  first  of  general  features  of  tho  data  input  and  operating 
procedure,  followed  by  details  of  specific  programmes  for  trajectory  and 
dynamic  analysis. 

Tho  format  of  data  input 

A  structural  format  has  been  used  for  all  programmes.  The  data  is 
grouped  into  five  tables  of  numbers;  a  table  of  atmosphere  properties;  a  table 
of  drag  properties  of  tho  body  or  of  moasured  accelerations  to  enable  tho 
trajectory  to  be  calculated,  a  list  of  variable  parameters,  a  list  of  constant 
parameters,  and  finally  the  experimental  data  to  le  analysed.  The  data  is 
headed  by  a  title  and  is  punched  on  five  hole  paper  tape. 

(a)  Atmosphere  table 

This  contains  measurements  of  temporature  and  pressure  made  at  the  time 
of  flight,  tho  units  being  °C  and  millibars.  Each  row  of  the  table  consists 
of: 

Height  (ft)  Pressure  (mbs)  Temperature  (°C) 
the  maximum  number  of  rows  is  28  and  the  table  is  terminated  by  the  warning 
character  (*),  an  asterisk.  The  programme  usos  a  log-linoar  interpolation  to 
obtain  density  and  Mach  number  as  a  function  of  height. 

(b)  Drag  table 

For  trajeotory  analysis  drag  coeffioiont  is  read  in  as  a  function  of 
Mach  number  (in  the  form:  Maah  number,  CA)  but  for  dynamic  analysis  height  is 
a  more  useful  variable  (vis:  height,  C^).  Wien  trajectory  analysis  makes  use 
of  in  flight  measurements  of  acceleration  this  table  takes  tho  form: 
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Time  (sec)  a  (g)  a  (g) 

a  y 

the  maximum  number  of  rows  is  32  and  the  tablo  is  terminated  by  an  asterisk  (*) 
warning  character. 

(c)  Variable  parameter  list 

In  the  present  programmes  up  to  12  parameters  may  be  variables.  The 
initial  estimates  of  the  parameters  are  read  as  a  table  along  with  the 
desired  perturbations  and  is  headed  by  the  initial  time,  thus: 

iv 


t 

0 


Xp  £p  * 

The  number  of  parameters  P  is  not  set,  instead  tho  table  i3  terminated  by  a 
warning  character  *.  If  desired,  parameters  may  bo  omitted  from  the  bottom 
of  the  list  when  they  are  assumed  tc  be  zero  by  the  programme  unless  also 
sot  in  the  list  of  constants. 

(d)  Constant  parameter  list 

This  is  a  list  of  all  constants  needed  by  the  dynamic  model  and  of 
cortain  numbers  which  control  the  programme.  Die  actual  list  depends  on  the 
particular  programme. 

(e)  Experimental  data  table 

IMs  table  is  headed  by  an  integer  q  which  is  one  loss  than  tho  number 
of  dimensions  of  the  experimental  data.  Hie  table  is  terminated  by  an 
asterisk  *.  If  desired  this  table  may  be  split  into  blocks,  headed  by  q  and 
terminated  by  *.  Tho  programme  then  fits  each  block  in  turn  and  uses  the 
value  of  the  parameters  obtained  for  a  block  as  the  initial  estimates  for  the 
subsequent  block.  In  the  general  three  dimensional  case  each  row  would  be 
wri tten  thus : 

\  ai  ^  Yi 

That  is  tho  timo  followed  by  the  components  of  the  observation  followed  by 
tho  weight.  The  value  of  q  in  this  case  would  be  2.  The  trajectory  programmes 
are  one  dimensional  so  that  q  is  zero  and  is  not  road  in.  The  dynamic  analysis 
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programmes  aro  either  one  or  two  dimensional  so  that  if  q  is  1  the  programme 
looks  for  and  (3^,  but  if  zero  only  for  a^.  In  those  programmes  is  not 
read  in  at  present  because  weights  are  assumed  to  be  unity. 

Operating  details 

The  programme  tape  is  first  road  into  the  computer  using  the  Mercury 
Autocode  Input  Scheme  or  the  Binary  Input.  The  machine  then  calls  and  roads 
tho  data  tape.  The  operation  of  the  programme  is  entirely  automatic  and  con¬ 
trolled  by  numbers  included  on  tho  data  tape.  Certain  automatic  controls, 
however,  can  be  over-ridden  using  the  hand-switches  as  will  bo  explained  below. 
For  normal  operation  only  hand-switch  4  is  sot  to  suppress  optional  printing, 
all  others  aro  cleared. 

(a)  Accuracy 

The  accuracy  to  which  tho  data  is  fitted  depends  on  the  number  of  itera¬ 
tions  mode  and  the  spoed  of  convergence.  Clearly  whon  all  the  changes  in  the 
parameters  are  all  zero  the  limit  of  accuracy  has  been  reached  but  this 
would  take  an  infinite  time.  Accordingly  it  is  arranged  that  whon  tho  bx^  are 
less  than  a  preset  level,  the  fitting  process  is  terminated  and  the  fitted 
motion  printod  out  (see  Fig.5)»  This  accuracy  levol  is  set  to  k  times  where 
aro  the  perturbations  being  made  in  the  parameters.  Those  perturbations  aro 
arbitrary  provided  they  are  small  onough  and  so  are  given  this  additional  use. 

After  each  iteration  the  new  value 3  of  the  parameters  are  printod  out 
along  with  the  variance  of  the  errors  and  tho  probable  error  in  the  poromotors. 
The  usor  is  thus  able  to  sec  how  fast  tho  solution  is  converging.  If  desired 
he  may  terminate  tho  fitting  process  at  ary  time  by  setting  hand-switch  1.  The 
last  fit  to  tho  data  is  then  printod  and  tho  programme  calls  in  the  next  block 
of  data. 

(b)  Rejection  of  data 

Some  moans  of  rejecting  data  that  is  obviously  wrong  is  necessary.  As 
mentioned  in  Section  6  this  is  done  by  rejecting  data  whon  tho  error  is  groator 
than  4  times  the  standard  deviation  of  all  errors.  On  the  first  fit,  however, 
this  lovel  a is  set  on  tho  data  tape  and  it  is  advisable  to  mako  it  as  largo 
as  possible,  say  about  ton  times  the  oxpected  errors  othorwise  data  is  rojoctod 
unnecessarily. 
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(o)  Interpolation  and  extrapolation  of  data 

Hie  programme  is  fitting  observations  at  discrete  time  intervals  but 
it  is  often  useful  to  have  the  fitted  motion  at  more  frequent  intervals. 

This  is  obtained  by  use  of  the  third  control  parameter  on  the  data  tape, 
the  maximum  step  length  h.  The  programme  will  print  out  the  fitted  motion 
at  this  interval  as  well  as  at  the  times  of  observations  so  that  interpolation 
is  easy.  If  h  is  greater  than  the  interval  between  two  observations  no 
intermediate  value  is  printed. 

To  extrapolate  before  the  first  observation,  it  is  only  necessary  to 

set  the  initial  time  to  the  appropriate  value.  The  motion  is  then  printed 

at  the  interval  h  up  to  the  time  of  the  first  observations.  To  extrapolate 

beyond  the  last  observation,  the  easiest  way  is  to  add  a  block  of  one 

observation  only,  (q  t  *)  will  do.  Hie  programme  will  then  extend  the 

fitted  motion  to  time  t  .  It  will  not  attempt  to  fit  a  block  of  less  than 

6 

P  observations,  it  will  merely  print  the  fitted  motion  using  the  latest 
values  of  parameters. 

(d)  Use  of  the  programmes  for  prediction 

The  facility  described  above  enables  the  programmes  to  be  used  for  pre¬ 
flight  predictions.  All  that  is  necessary  is  to  replace  the  experimental 

data  by  (0  t  0  *).  The  motion  between  t  and  t  will  then  be  printed  at 
e  o  e 

the  interval  h. 

Trajectory  analysis  programmes 

Trajectory  analysis  1623/lA.  This  programme  uses  a  table  of  Mach 
number  against  based  on  the  best  pre-flight  estimates.  The  parameter 
lists  are  as  follows: 

*0 

*CA 

Vo  £2  (ft/sec) 

ho  *  <ft> 

Vq  (ft/sec),  hQ  (ft),  0o  (deg),  s  (sq  ft),  m  (pounds),  Initial  Rejection 
error  e  (ft),  Accuracy  factor  k,  Max  step  length  h  (sec). 

The  data  fitted  consists  of  time  and  height,  the  weights  are  the  reciprocals 
of  the  probable  errors  of  each  ballistic  camera  height. 
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Trajectory  analysis  1623/lB.  This  programme  uses  measured  values  of 

acceleration  a  and  a  in  place  cf  the  drag  table  as  already  described.  The 
x  y 

parameter  lists  are  as  follows. 

t 

o 

A  ax  51  (ft/  sec2) 

V0  (fi/sec) 

ho  *  <ft> 

Vq  (ft/sec),  hQ  (ft),  6o  (deg),  s  (sq  ft),  m  (pounds),  mix  Cm(j/B,  CNa, 
e  (ft),  k,  h  (sec). 

Experimental  data,  which  is  height  from  ballistic  cameras,  is  the  same  as 
programme  1A.  Table  1  is  an  annotated  copy  of  data  as  read  into  the  computer 
and  Table  2  the  resulting  fitted  trajectory. 

Dynamic  analysis  programmes 

Four  programmes  are  available  at  present,  one  each  for  differential 
pressure  gauges  and  rate  gyro  analysis  and  two  for  accelerometer  analysis,  as 
follows. 

Dynamic  analysis  1623/1  -  For  differential  pressure  gauges  -  9  parameters 

1623/2  -  For  accelerometers  -  9  parameters 

1623/3  “  For  accelerometers  -  10  parameters 

1623/4  -  For  rate  gyros  -  9  parameters 

All  programmes  determine  the  parameters  C^,  C^,  P,  Rq,  Sq,  Yr  »  Ys 
and  two  zero  errors  for  programme  4  only,  and  two  trim  angles  in  the  case  of 
programmes  1-3.  At  present  the  latter  cannot  distinguish  between  zero  errors 
and  trim.  In  the  case  of  1623/3  the  extra  parameter  is  and  intended  to  be 
used  when  two  accelerometers  with  a  large  mutual  axial  displacement  are  being 
analysed.  The  parameter  lists  are  arranged  in  the  same  way  for  all  four  pro¬ 
grammes  and  are  as  follows:- 


°Mq  52 

P  (Radians  per  sec) 

C„  C,  Included  only  in  Prog.  1623/3 

R*  S5  (Rad) 

SQ  S6  (Rad) 
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Yr  (Rad) 

Ys  €g  (Rad) 

afc  (Rad) 

h  *10  (Rad) 


~o  rad/aeo~|  for  Prog.  1 6 23/4 

_e  rad/secJ 
y 


A  (slugs  ft  sq),  B  (slugs  ft  sq),  l  (ft),  s  (sq  ft),  m  (pounds),  C^, 

VQ  (ft/seo) ,  ho  (ft),  0o  (deg) 

e  (initial  rejection  error),  k  (accuracy  factor)  h  (sec,  max  step  length) 

x^/32‘2  y^/32*2  a^/32'2  ")  Accelerometer  coordinates  included  only  in 

x^32‘2  yg/32'2  z/32'2  J  Progs.  1623/2,  1623/3 

C^p  For  Prog.  1623/1  only.  The  differential  pressure  for  unit  dynamic 
pressure  at  one  radian  incidence  in  pounds  per  square  foot. 

Experimental  data  follows  the  standard  form  except  that  weights  are 
omitted,  either  one  instrument  or  an  ortnogonal  pair  may  be  analysed.  An 
illustration  of  the  use  of  Programme  1623/2  for  the  analysis  of  actual 
flight  data  is  given  in  Tables  3  and  4.  Table  3  is  an  annotated  version  of 
the  data  tape  and  Table  4  the  result. 
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SYMBOLS 

A  roll  M  of  I.  Also  for  acceleration 

B  pitch  M  of  I 

C  matrix  of  partial  derivatives 

C,  axial  force  coefficient 

A 

C.  lift  force  derivative 

La 

Cu  pitching  moment  derivative 

MCL 

C„  normal  force  derivative 

Na 

damping  in. pitch  derivative 
D  axial  force.  Also  matrix  of  residuals 
E  matrix  of  errors 

N  normal  force 

R  component  of  incidence.  Also  for  residual 

5  component  of  incidence.  Also  standard  deviation 

U  sum  of  squares 

V  velocity 

W  matrix  of  weights 

*k  ^  ParamSter  value 

a  constant,  also  acceleration 

b  constant 

f  partial  derivative 

g  partial  derivative  or  gravity 

l  reference  length 

m  mass 

n  number  of  observations 

P,p  spin,  number  of  parameters 

q,r  pitch  and  yaw  rate 

s  reference  area 

t  time 

weight  of  observation 
parameter 
observation 
¥  matrix 

a, (3  incidence  or  observation 

6  climb  angle 
perturbation 

p  air  density 

03  aerodynamic  frequency 
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Table  1 


TABLE  1 

TRAJECTORY  ANALYSIS  1633/ iB  DATA 


ATMOSPHERE 


H 

PRESS 

TEMP 

0 

998*1 

•3 

6000 

807*  X 

•4 

8000 

750*5 

XX 

*6750 

540.0 

-xo*5 

a  4000 

401*6 

“*4 

30000 

3x0*0 

“4* 

35000 

*45*5 

-5* 

4x000 

x8a*5 

-63 

43000 

167.5 

-64 

44000 

X59.3 

-6a 

5a  000 

106.9 

-68.5 

58000 

78.9 

-69*5 

65000 

55-7 

“63 

70000 

43*7 

“55 

75000 

34*5 

“53 

80000 

*7*4 

-48.5 

86000 

a  o*8 

“47 

9x000 

X6.7 

“47  * 

ACCELERATION  TABLE 

T(SECS)  AX 

AY 

8*5-745  -7*7 

5 

8*6.315  "xa*9 

7 

836.835  -17*6 

xo 

837*395  -34*0 

XX 

837*765  "33*1 

xa 

838*365  "44*o 

XI 

838*735  "53*0 

xo 

839*765  -65*5 

5 

830*3  55  -67*3 

a 

830*765  “65*9 

0 

831*365  -6a*x 

0 

831*765  -57*4 

0 

833*365  -53*6 

e 

833.765  -47.5 

0 

833*365  "41*8 

0 

833.765  "37** 

0 

834*365  -31.9 

0 

834*785  “*7*5 

0 

835*365  -33*6 

0 

VARIABLE  PARAMETERS 

8a6*x 

-30  5.0 

*S4S<>  5*® 

86500  15  • 


CONSTANT  PARAMETERS 

*545°  86500  88*4  >1185  xai  e*6 

5000  a.©  e*5 


1.8 
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BALLISTIC  CAMERA  DATA 


T (SECS) 

H(FT) 

S«RT  W 

8*6.5 

80467 

0.048 

8*6.6 

78959 

0.030 

836.7 

77473 

0.03* 

836.8 

75970 

0.036 

8*6.9 

74447 

0.033 

8*7.0 

73  947 

0.03* 

8*7. 1 

7*455 

o.o*8 

837.3 

69970 

o.o*6 

8*7.3 

68485 

0.0*4 

837.4 

6  7006 

0.0*  I 

837.5 

65550 

0.0*6 

8*7.6 

64089 

o.  0*0 

8*7.7 

6*648 

0.037 

8*7.8 

6 1*39 

o.ox6 

8*7.9 

59800 

0.0x9 

8s8.o 

58385 

0.0*0 

8*8,  t 

56983 

o.ox8 

8*8.3 

55603 

0.0x9 

838.3 

54*3° 

0.0*0 

838.4 

5*868 

0.0*3 

838.5 

5*5*5 

0.0*4 

8*8.6 

50x88 

0.0*4 

838.7 

48869 

0.0*3 

8*8.8 

47590 

0.03* 

8*8.9 

46309 

0.0*4 

8*9.0 

45046 

0.0*6 

8*9.1 

43801 

0.0*9 

8*9.3 

4*584 

0.0*4 

839.3  4*373  o.©»4 

839.4  40*37  0.018 

8*9.5  39088  o.ox* 

8*9.6  37968  0.009 


8*9.7 

36866 

0.008 

8*9.8 

35778 

o.  008 

839.9 

34699 

0.008 

830.0 

33644 

0.008 

830.x 

3*601 

0.009 

830.3 

31566 

O.OX  X 

830.3 

3055* 

o.ex6 

830.4 

*9553 

o.  o*6 

830.5 

38604 

0.030 

830.6 

3  76  8s 

0.0*8 

830.7 

*6778 

0.0*8 

830.8 

35894 

o.o*8 

830.5 

*5031 

0.030 

831.0 

34*06 

0.013 

831.x 

*339* 

0.0X4 

831.* 

33603 

0.0X4 

831.3 

3x8x3 

0.0x4 

83**4 

3x054 

0.0x3 

831.5 

*  o*  97 

0.0x3 

831.6 

*9589 

o.ox  X 

835.*  4300  x.o 
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Table  2 


TABLE  II 

TRAJECTORY  ANALYSIS  X633/1B  RESULT. 


-30 

5.0 

*5450 

5.0 

86500 

35.0 

0*46977 

49 

-63 

*•9 

Si# 

8.x 

8.6 

0.46975 

49 

-6a 

3.9 

*5466 

8.X 

866  5  x 

8.6 

TIME 

H 

8a6.x 

86651 

836.5 

836.6 

80503 

78978 

836.7 

77459 

836.8 

75945 

836.9 

74437 

8a  7.  « 

7*936 

837.*  7X44* 

837.3  69953 

8a 7.3  6*473 

837.4  670®! 

837.5  65538 

8a 7.6  64084 

837*7  63640 

837.8  6xao6 

837*9  597*4 

8a8.o  5*373 

838.1  56975 

8a 8.a  5559* 

838.3  54**© 

8a  8.4  5**64 

838.5  5*5*4 

8a 8.6  50*99 

838.7  48890 

8a8.8  4759* 

838.9  463*4 

8a9*o  45*67 

839.x  43*** 

839.3  43608 

839.3  4*407 

839.4  40**6 

839.5  39064 


SIQMA  N-P 
AX  ERROR  (FT/SEC. SQ) 
Vo  (FT/SEC) 


Ho  (FEET) 

V  AX  (G) 

DYN  P 

M  NO 

CA 

RES 

15466  -13.5 

7333 

15*6 

o.*54 

-36 

15380  -x6.o 

9409 

*5*5 

o.*53 

15338  — l6.8 

10070 

*5*5 

0*15* 

-19 

*5*7*  ”*7*6 

10789 

*5*5 

o.*47 

*4 

151x4  -18.3 

**5S* 

*5*5 

#•144 

*5 

15054  ”*9*5 

*»35* 

*5*4 

0.144 

to 

14989  -30.9 

13181 

*5*4 

••*45 

11 

*49*9  -»»*3 

*4054 

*5*3 

0.145 

*4 

14845  -33*6 

*4975 

*5*3 

0**44 

*7 

14767  -35.0 

16096 

*5*3 

•.14* 

is 

14683  -37.0 

17381 

*5*3 

0.143 

5 

*4593  ”»**9 

*8530 

*5*3 

0.143 

13 

*4497  “30-* 

19806 

*5*» 

0.143 

5 

• 

14395  “3*** 

3X117 

*5*» 

o.*43 

0 

14386 

-34.8 

33480 

*5** 

o.*43 

33 

*4*7* 

-37*0 

*3*93 

15.0 

©•*43 

16 

14048  -39** 

*5350 

15.0 

o.*43 

is 

0 

139*9  ”4**3 

36733 

14.8 

0*143 

8 

137*3 

“43*5 

38073 

*4*7 

o*X43 

13 

*3*39  ”45*5 

*943* 

*4*5 

0*143 

xo 

*34*9  ”47*» 

30803 

*4*3 

0*143 

4 

*3334  ”4**9 

331x8 

14.3 

0.X4* 

”9 

*3*74  “50.6 

333*7 

14.0 

0.14* 

—  11 

13008  -5*»3 

34496 

*3.7 

••141 

“3  1 

13837 

-53.8 

35649 

*3*5 

0.14© 

—0 

13663 

-55.* 

36775 

*3.3 

0.139 

“*5 

13483 

-56*4 

37866 

*3.* 

••138 

-3* 

13399  "57*7 

39006 

13*9 

0.137 

-37 

lain  "59*o 

40377 

13.7 

0.136 

“*4 

11919  ”60.3 

4*055 

**•5 

0.136 

”34 

11733  -6i.6 

4x667 

*3.3 

0.137 

x 

11533  -63.0 

4**07 

13. O 

O.X38 

*4 
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83  9*6 

37922 

113x7  -64.3 

42665 

11.7 

0.X40 

4* 

829.7 

36801 

1x108  -65.6 

43034 

11*4 

0.14* 

*5 

7* 

76 

829.8 

3570* 

*0895  -66.6 

433*o 

11.2 

0.143 

829.9 

34623 

10680  **66.9 

4343* 

10.9 

©•*43 

830.0 

335*7 

10464  -67.3 

43345 

xo.6 

0.144 

77 

830.1 

32532 

10246  -67.7 

43x82 

10.4 

0.146 

69 

830.2 

31518 

10028  *^68.0 

4*939 

XO.I 

0.147 

48 

830.3 

3<>5*7 

9808  -68.1 

426x8 

9*8 

0.148 

*4 

830.4 

*9557 

9590  -67.8 

42  xx6 

9»5 

0.150 

“4 

_ / 

830.5 

28610 

93  7*  “67.6 

4*4*5 

9»3 

0.151 

—0 

830.6 

27684 

9*54  “*7*3 

40687 

9.0 

0.154 

-a 

830.7 

26779 

8938  -67.0 

3990a 

8.8 

0.156 

-1 

830.8 

25897 

8723  -66.6 

39073 

8.5 

0.158 

“3 

830.9 

*5035 

85x0  -65.8 

382  12 

8.2 

0.160 

~4 

831.0 

24195 

8299  -65.0 

3  73*4 

8.0 

0.162 

x  x 

831.1 

23376 

8091  **64.3 

3*447 

7.8 

0.164 

*5 

831.2 

22  578 

7885  -63.5 

3555* 

7.6 

o.l66 

*  5 

831.3 

2 1800 

7682  -62.7 

34*33 

7*3 

0. 168 

*3 

831.4 

21042 

7481  -6i.8 

33694 

7.* 

0.170 

12 

831.5 

20304 

7284  -6o.8 

3*740 

6.9 

0.172 

“7 

831.6 

19585 

7090  ”59*9 

3*774 

6.7 

•  •*75 

4 

-6a 

7090 

19585 

831.6 

0.0 

0.0 

0.0 

19585 

7090  -59*9 

3*774 

*•7 

o.*75 

832.1 

16277 

6x64  -55.x 

26779 

5*8 

0.191 

832.6 

*34** 

53x6  -50.x 

21610 

4*9 

0.214 

833.1 

*0949 

455*  "44»* 

17005 

4*» 

0.242 

833.6 

8846 

3876  -39.6 

13x00 

3*5 

0.278 

834.1 

7061 

3*79  “34** 

9718 

2.9 

o.3*5 

834.6 

5555 

2761  -29.9 

7**7 

2.4 

0.383 

835.1 

4*89 

*3*4  “*5*8 

5228 

2  .o 

0.448 

835*» 

4062 

2232  -25.0 

4904 

2.0 

o.4*3 
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table  ill 

RE  ENTRY  DYNAMIC  ANALYSIS  16*3/*  DATA 

ATMOSPHERE  TABLE 
H  PRESS  TEMP 
87500  29*8  -49*0 

95Soo  *3*®  ”45*4 

205500  8,74  “44*7  * 

DRAQ  TABLE 
H  CA 

45000  o*x6a 
280000  0*155  * 

VARIABLE  PARAMETERS 
787*680 
-0*30  0*01 

-a •©  x*o 
33*0  0.2 

0*06  o*oox 
0*05  0*002 
-o*6  0*03 

3*7  0*03 

O  0*002 
O  0*002  * 

CONSTANT  PARAMETERS 

0*520  3*365  2*333  2*4  224*5  2*85 

24335*7  96394*4  84*0  3*©  3.0  0*05 

ACCELEROMETER  COORDS (FT/Q) 

-0*00985  0*00560  -0*00268 

0*00300  0*00543  -0*00328 


ACCELEROMETER 

DATA 

time 

Y  ACC 

Z  ACC 

X 

787.700 

8*24 

-5*34 

787*730 

4*64 

-9.45 

787.740 

-X.79 

-8*32 

787*760 

-4*3* 

-2. 52 

787*780 

0*50 

5*04 

787*800 

7*26 

6*04 

787*830 

8*38 

0*76 

787*840 

2*53 

-3*93 

787*860 

-6.82 

-0*78 

787*880 

-20*07 

5.90 

787*890 

-0*30 

9.05 

787*900 

-5.96 

20*38 

787*960 

-9.83 

-4*44 

787*970 

-23*33 

-4*28 

788*000 

-23*75 

3*3  0 

788*020 

-20*37 

3*30 

788*0*0 

-6.85 

*  *44 

788*040 

-0*39 

-4.5* 

788*060 

-3.56 

—  S3  *08 

788*080 

-9*70 

-22*63 

788*200 

-9*80 

-4.5* 

788*210 

0*72 

O.73 

788*240 

8.58 

-2*48 

788*260 

8.45 

-7*27 

788*280 

2*09 

-6*40 
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TABLE  IV 

RE  ENTRY  DYNAMICS  ANALYSIS  16*2/*  RESULT, 


”9,30000 

0.0X00 

”3.00000 

I. 0000 

*3,00000 

0. xooo 
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*7 
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X.4500 
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0.07555 
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0.047* 
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0.00*0 
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4* 
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FIG.  I  DEFINITION  OF  NOTATION  &  FORCES 

ON  MISSILE 
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Fig. 3 


Fig.3  A  Black  Knight  instrumentation  assembly 
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FIG.  5  FLOW  DIAGRAM  OF  AUTOMATIC 
ANALYSIS  PROGRAMME 
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Fig. 7 


FIG. 7  FITTED  CURVE  TO  DATA  FROM  PAIR 
OF  ACCELEROMETERS  (2R~S) 
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FIG.  8  FITTED  CURVE  TO  SINGLE  RATE  GYRO  RECORD 
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FIG. 9  VARIATION  OF  PARAMETER  VALUE  WITH 
ITERATION  No  FOR  AN  ACCELEROMETER  ANALYSIS 
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